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ABSTRACT 


A three-dimensional numerical model of the Jet Engine Hush House located at Na- 
val Air Station Jacksonville. Florida was developed using the PLIOENICS computer 
meee ihe PFLOENICS code was used to determine the steady state aerothermal char- 
acteristics in the Flush Flouse produced by a Navy F-4 (Phantom IT) J-79-GE-$ gas 
iMeeire engme with afterburner. The PHOENICS generated pressure, velocity. and 
temperature fields of the test cell facility based on the k-e turbulence model are com- 
pared with results generated by a model incorporating a temperature dependent 
kinematic turbulent viscosity. In addition, for the k-e turbulence model. the effect of 
augmenter tube length on svtem behavior was investigated. Lessons learned and rec- 
onimendations in obtaining a converged solution are included to aid in further applica- 


tions of the PHOENICS code. 
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I. INFRODUCTION 


Todavs sophisticated aircraft engines, both new and overhauled. require full opera- 
tional testing prior to being placed in service. Jet engine test facilities are utilized to 
safelv and effectively monitor aircraft and jet engine performance characteristics. In 
addition to providing an environment for acquiring test data, the test cells supply clean, 
distortion-free airflow to the engine inlet, decrease noise pollution by decreasing the 
kinetic energy of the jet, and cool the exhaust prior to exiting to the atmosphere [Ref. 
Pa. ||. 


Miitereare five types of tests facilities currently used bv the NAVY. 


1. The Run-up Pad--where the aircraft 1s tied down in an unprotected environment 
to conduct power checks. 


2 


Sound Suppressors--equipment packages that reduce the noise of aircraft on the 
Run-Up Pad. 


Open Test Stand (OTS)--outdoor engine test stand and portable sensing equip- 
ment. 


2 


4. Hush House--a total aircraft enclosure capable of around the clock testing. 


can 


Fully Enclosed Test Cell(TC)--an enclosed facility for complex testing of engines 
out of the airframe. 


The run-up pad. sound suppressor. and the open test stand facilities although relatively 
inexpensive are highly inefficient and do not support military readiness needs. Noise and 
air pollution factors determine the allowable locations of the facilities and the environ- 
niental conditions control the time of operation of some of the facilities. The hush house 
and the fully enclosed test cell can operate independent from environmental conditions, 
and since noise is reduced, around the clock testing is allowed. [Ref. 2: p. 2] 

As of March 1984, 75 Navy and Marine activities have 76 fully enclosed test cells 
and five hush houses. Only 60 out of the 76 test cells are fully operational. Most of the 


test cells are over 27 vears old. Modern jet engines operate at much higher temperatures, 


mass flows, and velocities than thev did 20 vears ago when most of the test cells were 
designed and built. Due to the age of the facilities and high gas temperatures and water 
impingement. cracking and spalling of the test cell concrete walls continues to occur. 
In addition to the structural degradation of the test facilities many of the test cells are 
obsolete in that they are unable to support the increased size and wide variety of the 
engines requiring testing. [Ref 2: pp. 3-4] 

Noise, air, and water pollution control at test facilities 1s under increasing scrutiny. 
Prior to 1976, military test cells were exempt from state and local regulations pursuant 
to section 233 of the Clean Air Act. The Federal Court decided in 198], that jesttengam. 
test cells were to adhere to the requirements of the Clean Air Act. The absences 
degradation of exsisting pollution control systems has increased the risk of physical in- 
jury to the test cell facility personnel as well as damage risk to test cells, equipment, and 
aircraft. In 198i, the Navy Environmental Health Center cited hearnng losssamdiiesame 
atory damage for $114 million with an anucipated claim growth of 99% each vear. 
Equipment and aireraft are subjected to acoustic fatigue which could possibly shorten 
Operating life Rel sn) 

It is the Navy's goal to provide standardized facilities for testing aircraft and aircraft 
jet engines that will incorporate the latest technology and which will adhere to all local. 
state, and federal regulations concerning environinental safety. The current jet engine 
test cell facilities were designed empirically. based upon the designer’s experience and 
occasional aerothermal data obtained from existing test cells. The test cell facilities in 
general performed well for the type engines that they were designed for. Due to the in- 
crease in complexity and power of modern jet engines and aircraft the test cell facility 
design procedures of the past are inadequate. The design of a modern jet engine test 
facilitv must incorporate the effects of multiple engines producing nonconcentric flows, 


turbulence, supersonic velocities. extreme temperature gradients, and awkward 


gseometries. In addition. the new test cells must be able to accomodate a wide variety 
of engines and aircraft each with its own characteristic exhaust temperature. exit 
velocity, and mass flow rate [Ref 1]: p. 1]. 

The design of modern test cell facilities requires a more sophisticated design proce- 
dure. A numerical simulation of the test cell can provide the design engineer with criucal 
performance data from existing test cells to aid in the design of new test cell facilities. 
The numerical model provides a cost and time efficient means of obtaining the per- 
formance data of existing test cell facility designs, and this information can be incorpo- 
rated in the design of new test facilities. 

The PHOENICS Code ts used to numerically model the aerothermal characteristics 
in the Hush House Facility. located at Naval Air Station Jacksonville. Flonda, due to 


the operation of a Navy J-79-GE-S gas turbine engine. 


II. PROBLEM DESCRIPTION 


A. HUSH HOUSE CHARACTERIZATION 

In 1980 the Air Force purchased 79 standardized demountable Hush Houses with 
test cell capabilities. The large scale purchase, standardized design, and simple 
contruction resulted in a relatively low cost of $2.2 million per Hush House. The Hush 
House located at Naval Air Station Jacksonville, Florida is of the same design as the 
Air Force Hush House except for the additional AF FF Firefighting svstem tstalled bv 
the Navy. Neither the Air Force’s nor the Navv’s Hush Flouse design incorporates air 
pollution control] systems. Although the cost of the initial large scale purchase bv the 
Air Force was low, the overall long term cost was significantiv Iigher due to tam 
maintenance repair costs. It appears that the purchase was made in haste and not fully 
researched. The Navy could benefit from a similar large scale purchase providing that 
a high performance low cost maintenance svstem can be designed. [Ref. 2: P. 1] 

The basic operation of the Hush House consists of four major sections: intakes, test 
bay, augmenter tube, and exhaust stack as shown in Figure | on page 6 and Figure 2 
on page 7. The primary intakes provide air for the gas turbine engine intake. The sec- 
ondary intakes provide air to mix with the engine exhaust. The mixing of the low tem- 
perature low velocity secondary air with the high temperature high velocity exhaust 
results in cooling the hot exhaust gas and reducing the exhaust gas velocity. The test 
bay is a large unobstructed area where the entire aircraft or an engine mounted in a test 
stand 1s positioned for testing. The large test bav allows intake air to smooth out into 
a uniform profile prior to the engine intake thereby approximating true operating con- 
ditions. The augmenter tube is lined with acoustic absorption material to aid in acoustic 


energy dissipation from the engine and jet flow. The augmenter tube aids in the mixing 


of the engine exhaust and secondary intake cooling air. The exhaust stack changes the 
flow direction from horizontal to vertical as the exhaust gas flow exits to the atmos- 
mmere. |Kef. |: p. 2| 

The overall dimensions of the test bay are 28.7 meters (94 feet) wide bv 25.5 meters 
($3.7 feet) long by 8.6 meters (28.2 feet) high. The 24.4 meter (76 feet) long oval shaped 
augmenter tube is centered at the rear end of the test bay with the floor of the tube at 
the same level as the test bay floor. The oval augmenter tube cross section is roughly 
4.9 meters (16 feet) wide by 4.0 meters (13.2 feet) high. The dimensions are constant 
throughout the tube length. The exhaust stack at the end of the augmenter tube 1s 
simplv a 45 degree ramp ending just short of the top of the augmenter tube. The exhaust 
Meee 1s a $.2 meter (27 feet) by 8.2 meter (27 feet) horizontal plane. The aircraft is 
positioned in the test bay with the engine outlet approximately 4.6 meters (15 feet) from 
the augmenter tube entrance (see Figure | on page 6 and Figure 2 on page 7). General 


Operating procedures for conducting engine tests are described bv Nicolaus (1988) [Ref. 
3: pp. 5-6]. 


B. ISSUES CONCERNING THE HUSH HOUSE 

The Air Force type Hush House in general is performing satisfactorly. Some struc- 
tural and operational problems do exist in the Hush House. The stainless steel corru- 
gated sheets that overlay the acoustic absorption material is subject to cracking along 
the welds. During long periods of idle engine operation a portion of the exhaust gas flow 
is deflected away from the augmenter tube entrance. The exhaust gas rises to the top 
of the test bay activating the high heat sensors [Ref. 3: p. 7]. Vibration of the entire 
Hush House structure is present during engine testing. It is not exactly clear what effect 
the vibration has on the structural integrity of the Hush Houses. Some of the Hush 
Houses require inspection and retightening of the bolts that fasten the acoustic absorp- 


tion material to the walls of the augmenter tube [Ref. 2: p. 4]. 
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Figure 1. Actual Hush House (Top Veiw Looking Down) 
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Figure 2. Actual Hush House (Side View) 


Hl, THE PHOENICS CODE 


A. GENERAL DESCRIPTION 

PHOENICS 1s a commercially available general purpose computer program copy 
nighted bv CHAM (Concentration Heat and Momentum) Limited which was founded 
in 1972 bv Prof. D. Brian Spalding of Imperial College in London, England. 
PHOENICS, an acronym for Parabolic Hyperbolic Or Elliptic Numerical Integration 
Code Series, simulates fluid flow, heat transfer, and chemical reactions.[Ref. 4: p. 1.1] 

The mathematical basis of PHOENICS 1s the iterative solution of descretized con- 
servation equations of mass, chemical species, momentum, energy. turbulence quantities, 
and other related fluid properties. The general form of the conservation equation (before 


descretization) for which PHOENICS solves can be expressed in the form |Ref. 85 paeaae 


“~ 


7 (p@) + dir(pl'd = Pgrad,) = Sy (3.1) 
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where: 


@ = dependent variable 


a 
oO 
| 


exchange coefhicient for @ 


Y 
© 


= source of @ per Unit vollime 


p = densitv 


= 


VelOCciO VEELor 


ro 


b= 


example, in the energy conservation equation the dependent variable @ would be 
enthalpy and the exchange coefficient T., would be the effective thermal conductivty di- 
vided by the constant pressure specific heat. The partial differential equations(p.d.e.) 
are coupled and nonlinear. The coupled, nonlinear conservation p.d.e’s can only be 


solved (except for some Very simple cases ) by approximate numerical methods. Typical 


methods for their approximate numerical solution includes the finite difference and finite 
element methods. In the present case, the control volume method has features of both 
the finite difference and finite element methods. All these methods transform a partial 


differential equation into system of algebraic equations as follows: 


ao» = AvOy+ ADs a arQ 5 =e Ay? wor Ay@ 1 =e arPy, Sie azirt+ b (ee) 


where: 
a = influence coefficients 


b 


lI 


a representation of the source of @ for the cell 


ieee subscripts N, S$. E,W, HH, and L denote the North, South, East.West, High, and 
Low neighboring cell faces respectively as shown in Figure 3 on page 13. The subscript 
P denotes the grid node at a particular location and the T subscript denotes the same 
grid node at an earlier tume. The system domain 1s partitioned into a grid network of 
cells. This partition discretizes the continous domain into .V. cells. Then each p.d.e. is 


transformed into an algebraic equation for each cell as follows: 


b _ AvOy. + Ae ar@rt AnD yt AryO7, + Q,0, + aror+ Ss | 
= : 


Ax + a as Gp ap Gis Crain Gia art Ap 


oe) 
(42 
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Thus for example, the energy p.d.e. 1s transformed into .\; algebraic equations which are 
then solved for the value of enthalpy ( or temperature ) at the center of each of the \; 
cells. It should be noted that these algebraic equations are also coupled, and therefore 
cannot be solved without solving the discretized algebraic equations for all of the de- 
pendent variables. The algebraic equations are numerous and dependent on similar 
equations from neighboring cells. The number of algebraic equations is equa] to the 
number of dependent variables times the number of cells, VN... The Hush House utilizes 


4800 cells. Solving for seven dependent variables in each cell equates to the solution of 


33600 equations. PHOENICS uses the SIMPLE algorithm requiring a guess and cor- 
rection scheme of which the objective is to reduce the imbalance between the left and 
right hand sides of the discretized equation to a negligible magnitude [Ref. 5: p. 2.5]. 
The specific methods of reducing the rmbalances (referred to as residuals) and thereby 
obtaining a converged solution 1s discussed in Chapter four. 

The PHOENICS computer code consists of two main programs EARTH and SAT- 
ELLITE. Figure 4 on page 14 shows the relationships between EARTH, SAT EUEDDE 
and their respective subprograms and hardware. The EARTH program 1s the largest of 
the two and it contams the coding for solving the conservation and related phvsical law 
equations. SATELLITE-EARTHI communication ts uni-directional, information can go 
onlv from SATELLITE to EARTH. Therefore a property which itself is a funeno@imen 
the solution, which is achieved in EARTH. cannot be accommodated by coding onlv in 
SATELLITE. For accommodating solution dependent properties. one must use 
GROUND. The subroutine GRENI is an exemplary GROUND code that)previds: 
common options for general material properties such as density and viscosity. Desired 
options which are not in GRENXI must be programmed in GROUND. 

The SATELLITE program transposes input data expressed in PHOENICS input 
language (PIL) into commands that will actvate EARTH. A QI file (quick-input file. 
number 1), where the flow problem ts defined, 1s a formatted input file divided into 24 
groups: 

GROUP 1 = Run utle and other preliminaries 
GROUP 2 = Transcience: time-step specification 
GROUP 3 = N-direction grid specification 

GROUP 4 = Y-direction grid specification 

GROUP 5 = Z-direction grid specification 

GROUP 6 = Body-fitted coordinates or grid distortion 


GROUP 7 = Variables stored, solved and named 


GROUP &§ = Terms (in differential equations) and devices 


GROUP 9 = Properties of the medium (or media) 
GROUP 10 = Iterphase-transfer prosses and properties 
GROUP 11 = Iniualization of variable or porosity fields 
GROUP 12 = Convection and diffusion adjustments 
GROUP 13 = Boundary conditions and special sources 
GROUP 14 = downstream pressure for PARAB=T 
GROUP 15 = Termination of sweeps 

GROUP 16 = Termination of iterations 

GISOUP |7 = Underrelaxation devices 

GROUP 18 = Limits on variables or increments to them 
wc OUP 19 = Data communication by SATELLITE to GROUND 
GROUP 20 = Prehminarv printout 

GROLP 21 = Print out of variables 

GROLP 22 = Spot-value printout 

GROUP 25 = Field printout and plot control 

GROUP 24 = Dumps for restarts 


A complete explanation of each group’s data structure and a PHOENICS Input Lan- 
guage glossary can be found in the PHOENICS Beginners Guide and Users Manual (see 


Feorerence 2). 


B. MODELING THE PROBLEM 

The main scope of the thesis is a parametric study of the Naval Air Station, 
Jacksonville. Florida, Hush House originally modeled by Nicolaus in 1988 [Ref. 3]. The 
original model of the Hush House used the cartesian coordinate grid system of X, Y. and 
Z coordinates. The original three dimensional rectangular grid consisted of 40 cells in 
the Z direction (NZ=40, low to high), 12 cells in the Y direction (NY= 12, south to 
north), and 10 cells in the X direction (NX = 10, west to east). The Hush House model 
takes advantage of the Y-Z plane of symmetry of the actual Hush House and therefore 


Om walt or the Hlush House is modeled. Dividine the structure in half reduces the 


1] 


nuinber of equations to be solved bv half thereby reducing computer CPU time. The test 
bay is defined by cells 1-19 in the Z direction, cells I-12 in the Y direction, andwaememe] 
in the X direction. The jet engine is defined by cells 12-14 m the Z direction. celigmaas 
in the ¥ direction, and cell 1 in the X direction. Ihe augmenter tube 1s deiimedsiueenls 
1$-27 in the Z direction, cells 1-8 in the Y direction, and cells 1-5 in the X direction. The 
exhaust stack is defined bv cells 28-39 in the Z direction , cells 1-8 in the Y direction, and 
cells 1-5 in the X direction. The 45 degree ramp at the end of the exhaust 1s modeled 
with a step svstem. Figure 5 on page 15 and Figure 6 on page 16 graphically display 
the \-Z and the Y-Z planes of the structures |Rel. oppo! 

The graphical results of the temperature and velocity profiles at volume slab IX= 1 
for IY =4 from the original model shown in Figure 7 on page 17 and Figure § on page 
18 indicate that the exhaust gas temperature and velocity were reduced bv roughlv 30 
percent in the first 20 percent of the augmenter tube [Ref. 3: pp. 35-356), Thevquesiaam 
then arose as to what would be the effect of a shorter auecmenter tube) fies 
parametric studv investigates the effects of augmenter tube length on the cverall per- 
formance of the Hush House. Recall that the purpose of the augmenter tube 1s to re- 
duce the exhaust gas temperature and velocity to an environmentelv acceptable level 
prior to releasing the gas to the atmosphere. The performance of the Hush House 1s a 


function of the following parameters: 


]. Condition of the flow field prior to entering the engine inlet 


tJ 


Magnitude of the mass flow rate of the secondarv cooling air 


3. Reduction of the gas temperature and velocity through the augmenter tube 


The original augmenter tube length of 21.0 meters was reduced 10% to a length of 
18.9 meters. The parametric study on the effects of augmenter tube length on svstem 


behavior used the same k- e turbulence model as the original model. If the performance 





Figure 3. Cell Faces 


of the IJush House was not significantly reduced by a shorter augmenter tube then the 
cost and time of construction, as well as the building size could be reduced. 

Prior to conducting further geometrically based parametric studies the question 
arose of the effects that different turbulence models had on the aerothermal character- 
istics generated by PHOENICS. The concern is to know how dependent the pressure, 
velocity, and temperature fields were on turbulence modeling. The GRENI subroutine 


@ierseviye uscr of PHOENICS four turbulence models: [Kef. 4: p. 3] 
1. Algebraic Model 


ty 


Mixing Length Model 
k-Mixing Length Model 
. k- € Model 


bow 


Gas flow through the IIush Tlouse is sinular to flow through a three dimensional 
duct which is a wall boundary problem between two parallel Y-Z planes and two parallel 


X-Z planes. The Hush Flouse model is actually much more complicated then just stated 
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but in attempting to select an existing turbulence model an assumption of the geometry 
configuration is required. The algebraic and mixing length models as coded in GRENI 
are applicable to flow over a single flat plate and therefore are not suitable for this 


model. In the k-mixing length model turbulent viscosity is calculated according to: 
0.5 | 
= Con G (3.4) 


where k, the turbulence kinetic energy is a dependent variable solved for by PHOENICS. 
(yas a constant set in GREAI and @, is the mixing length. GREN1 has seven options 
for computing the mixing length /,. The first three options. are applicable to flow over 
a single flat plate and are not suitable for the Hush House model. The fourth option 
selects the mixing length scale suitable for free jets. The Hush House has a wall bounded 
jet. The fifth option is suitable for free jets, but requires the setting of nine constants 
which can only be determined through experimentation of the specific flow svstem. The 
unavailability of these constants ruled this option out. The remaining two options select 
Nikuradse’s mixing length seale which requires a cylindrical polar coordinate grid system 
for flow through a cvlindrical duct. The k- ¢ turbulence model ts the onlv PHOENICS 
model that is suitable for the Hush House without drastic changes to grid geometry. 
Pi@svever, since the k- ¢ turbulence model was used originally and the other GRENI 
turbulence options are unsuitable a turbulence model had to be developed and inserted 
into GROUND. A simple temperature dependent turbulent viscosity model was chosen 
because of the complex geometry involed. The turbulent viscosity term ENUT was set 


equal to a multiplying factor times the kinematic laminar viscositv ENUL. 
fone | er x EN (3.3) 


ENUL (commonly referred to as v,) is a function of the dvnamic viscosity uw and the 


MENSILy 9p . 
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ane Lu oe 
bales = a (3.6) 
The dynamic viscosity of air is temperature dependent over large temperature ranges. 


u(295 K) = 18.9 x 10° HS aun 


Hr 


(200° K) = 62.27 x 1976 AS (3.8) 


aa’ 


A fourth order polynonual equation for u (T) was developed by Andrews and Biblarz in 


1981 [Ref. 6: p.C-1]. 


~~" @ 


u(T) = (xl +x2x T+x3x To +xdx To 4x5x T) x 10x 107 (3.9) 


where: 
x1 = 4.8856745 
\2 Sees = 10 
x3 = —2.4261775 x 10> 
x4 = 7.9306 x 10-9 


cm 
lI 


X NUS 9S Sal al 


PHOENICS solves for the temperature and density of cach cell therefore Naas 


temperature-density dependent and is not a constant value for each cel] in the grid for 


anv given multiplving factor, 


[EMOn Pe Nee = (3.10) 





The code line ENUT=GRND in group 9 of the QI file directs°SA VTEERIVE tomme 


GROUND subroutine to evaluate ENUT for each cell. The code for evaluatine ENG 


is imbedded in group 9, section 6 of the GROUND subroutine enclosed in Appendix 
cS. 

The multiplying factor was arbitrarily selected as 200 initially and would be in- 
creased by 200 after convergence of each solution. This procedure would produce sets 


of pressure, velocity. and temperature fields that are based on different kinamatic tur- 


- 


Lye 


bulent viscosity values ranging from .0001 to 10 re It was anticipated that for low 
values of kinematic turbulent viscosity the mixing effect of turbulence would be low and 
therefore the flow field would be highly streamlined. As the turbulent viscosity and 
mixing effect increased it was anticipated that the momentum transport between the 


exhaust gas and entrained cooling air would increase, resulting in lower velocities and 


temperatures exiting the exhaust stack. 


IV. SOLUTION CONVERGENCE 


A. SOLUTION METHODS 
PHOENICS provides three methods for solving the NA x NY x NZ nUiieomeen 


equations for each dependent variable: 


1. Slabwise Method 
2) Miiole Field Slethod 


3. Point-Bv-Point Method 


A slab 1s an array of cells with the same value of low-to-high in the Z coordinate. 
In the Hush House model there are 40 slabs. The Slabwise method of solution solves for 
all the dependent variables in the lowest slab and then sequentially continues to the 
highest slab. The sequence of solving equations from low to high slabs is referred to as 
a sweep. The solution of the discretized equation for any dependent variable requires 
values from each neighboring cell. In the Slabwise solution the equations for values in 
a slab require values from the slabs below and above it. Adjustments to the values in 
the next higher slab will invalidate the values just solved for in the lower slab, therefore 
manv sweeps are required until all the discretized equations balance and no further ad- 
justmenits to nershboring cells are required Kel aap eo) 

The Whole Field Method of solution simultaneously solves the NX x NY x NZ 
number of equations for each dependent variable in one sweep thereby reducing the 
number of sweeps compared with the Slabwise Method bv a factor of NZ. The Whole 
Field Method requires much more computer storage space then the Slabwise Method 
and should not be used if storage space 1s hmited. [Ref. 7: p. 2.5] 

The Point-By-Point Method of solution does not employ PHOENICS simultaneous 


linear equation solver for the evaluation of the dependent variables. The node values 


to 
tl 


of each variable are updated by simple arithmetic substitution. This simple substitution 
reduces the amount of computational work but requires many cycles to reach a final 


solution. 


B. METHODS OF DETERMINING CONVERGENCE 

A converged solution 1s one in which the imbalance between the left and right hand 
sides of every discretized equation is reduced to an acceptable small value. Due to the 
iterative procedure utilized by PHOENICS in solving the discretized equations. conver- 
gence is not guaranteed. User intervention is usually required to obtain a converged 
solution. 

There are four means available to the user of PHIOENICS to evaluate whether a 


solution is converging or diverging: 


1. Tabulation of residuals vs sweep 
2. Conservation of mass flow 
3. Conservation of energy 


4. Comparison of successive sweeps 


Residuals calculated for each dependent variable are the sum of the absolute values of 
the discretized equation imbalance from each cell. The decrease in residuals with each 
successive sweep 1s the primary means of determining convergence. The output data file 
generated by PHOENICS provides the values of residuals vs sweeps for each dependent 
Variable in tabular and graphical form. Figure 9 on page 28 show graphical examples 
indicating monotonic convergence. Monotonic convergence is characterized by sequen- 
tial decrease of residuals vs sweeps. Monotonic convergence will usually, but not al- 
ways, vield a final solution in the shortest amount of time. Monotonic convergence 1s 
not alwavs possible or may be very difficult to obtain. Residuals mav have a damped 
oscillatory decay as shown in Figure 10 on page 29. The computational time required 


to reach a converged solution with the residuals converging this way may be 


considerably greater than 1f monotonic convergence could be obtained. [lowever a 
rapidly decreasing damped oscillatory decay can converge faster than o slowly decreasing 
monotonic decay. Divergence of a solution is determined by the residuals increasing 
with the number of sweeps as shown in Figure 1] on page 30. A diverging solution if 
unaltered will most likelv result in “self destruction” of the simulation. 

There are no cut and dry methods for obtaining and maintaining a converged sol- 
ution. The general guidelines for obtaining a converged solution are discussed in part 2 
section 8 of Reference 7. In choosing a viable set of convergence control options from 
the large number of possibilities, the PHOENICS programmer must use common sense, 
experience, and trial and error. 


Three convergence control options were utilized in the Hush House model: 


1. Under relaxation factor 
2. Ternunation criteria for mner iterations 


3. Incrementation of source terms 


The under relaxation factor controls the rate of change of a dependent variable from 
sweep to sweep. Due to the iterative procedure utilized by PHOENICS, impropemean 
ting of the under relaxation factor of a dependent variable can result in a diverging sol- 
ution. Under relaxation factors range from just greater than zero to one. Small values 
of the under relaxation factor require more iterations for convergence to be obtained and 
thus more computer time to reach a final solution. An optimum under relaxation factor 
1s one that minimizes computer time and also prevents divergence. The determination 
of an optimum under relaxation factor, as in the choice of convergence control options, 
requires experience and trial and error. The under relaxation factor is activated by the 
use of the RELAX command in group 17 of the QI! file. 

The second convergence option, termination criteria for inner iterations, controls the 


number of iterations performed on a dependent variable for each sweep. A larger 


number of iterations per sweep will vield a more accurate value for that variable and thus 
result in fewer sweeps but the computing time per sweep is increased. <A balance 
between the number of inner iterations per sweep and the number of sweeps is required. 
It would be a waste of computing time to iterate on a variable to a degree of accuracy 
exceeding the required accuracy of the output. The number of inner iterations is 
controled by the command LITER in group 16 of the Q} file. 

The incrementation of source terms is a method used to slowly apply the excitation 
to the flow svstem thereby reducing the jump required in solving the conservation 
equations. The total excitation to the system is obtained by applving the source, such 
as mass flow rate or temperature, in small increments. The magnitude of the excitation 
1s increased after run of sweeps has converged based on the previous excitation magni- 
tude. This procedure is in effect a transient to steady state flow condition solution. 

For the first parametric study involving the effects of augmenter tube length on 
system behavior, a converged solution with a ten percent reduction in augmenter tube 
leet was desired. Having little experience at this time with PHOENICS, the same 
under relaxation factors and number of inner iterations were used as in the original Hush 
mise model ©] file developed by Nicolaus. The variable FACT 1s used to define the 
magnitude of the FALSDT under relaxation factor. A value of FACT = 0.75 was ap- 
@ieanto ihe dependent variables U1, VI, WI, KE, and EP. A value of 2.5 x FACT was 
used for the dependent variable [1]. The variable LINRLX 1s used to define the linear 
under relaxation factor. A value of LINRLX = 0.3 was used for the variable P1. The 
default value of LITER equal to 20 iterations per sweep was used for all dependent 
variables. The DFO09 file from the original Hlush House converged solution was re- 
started. To avoid a jump or shock to the system the original augmenter tube length was 
initially reduced by 5 percent with 2500 sweeps required for a converged solution. The 


augmenter tube length was shortened an additional 5 percent of the original length with 


3500 sweeps required for a converged solution. Each of the 32 runs, consisting of either 
300 or SO sweeps, exhibited monotonic convergence. 

For the second parametric study involving the effects of turbulent viscositv on svs- 
tem behavior a trial and error procedure for obtaining convergence was required. The 
multiplying factor MF of equation (3.10) was set equal to 200. The DF09 file from the 
original Hush House model was used for a restart with the exceptions of the dependent 
variables KE and EP which were not required for the new temperature dependent tur- 
bulence model] . Still having limited experience with PHOENICS, the same under relax- 
ation factors and number of inner iterations were used as in the original Hush House 
model Q] file. The above procedure resulted in divergence of the solution. The number 
of iterations, LITER, for the pressure variable Pl was increased from the default value 
of 20 to 40 because in the SIMPLE algorithim the pressure is solved for first from the 
continuity equation and then the velocitv components are solved from the momentum 
equations Which have a pressure gradient term. A more accurate value of pressure ob- 
tained by more iterations will vield more accurate velocity components without having 
to increase the number of iterations performed on them. The variable LINRLX was 
reduced from 0,3 to 0.1 for pressure variable Pl. The under relaxation factor) Fae 
was reduced from 0.75 to 0.25 for the remaining variables L1, V1, W1, HI. Dheseaminmes 
changes resulted in an oseillation of residuals with an overall mean effect of neither in- 
creasing nor decreasing. [It was evident that after nearly a thousand sweeps that the 
solution was not going to converge within an acceptable time frame. The next attempt 
at obtaining a converged solution involved incrementing the temperature source from 
the jet engine exhaust. The exhaust temp was set equal to ambient air at 295 ° Kk and 
the problem was started from scratch and not a restart from the original model. The 
niass flow rate was maintained at 38.49 kg s. The first few runs, consisting of $0 sweeps 


each, exhibited monotonic convergence. FACT was then inereased from 0.25 to 1.0 to 


minimize computer CPL ume. Monotonic convergence was present during the 1680 
sweeps required for a final solution with engine exhaust temperature at 295 ° kK. Engine 
exhaust temperature was increased to 500 ° K and a converged solution was obtained in 
5280 sweeps. With the engine exhaust temperature incremented to 700 ° Kk, an addi- 
tional 3100 sweeps Were required for a converged solution. With the engine exhaust 
Bemiperature set to 900 © Kk. the residuals oscillated severly, vet the mean values were 
neither converging nor diverging over successive runs. FACT was incremently reduced 
to 0.1 and LINRLA remained at 0.1. After nearly 5000 sweeps, oscillation of the resi- 
duals continued with the same mean effect as exhibited by the residuals of the first few 
runs. This dilemma of non convergence non divergence 1s a prime example of not 
choosing the correct convergence contro! options. It should be fullv understood that the 
solution of the flow field with engine exhaust temperature at 900 ° K 1s feasible but it 
requires a different approach than that used for the lower temperatures. Due to the 
shortage of time remaining for this turbulence study and the anticipated convergence 
problems at higher temperatures. it was decided to conduct the parametric study of tur- 
bulence using the converged solution with engine exhaust temperature at 700 ° k. 
Three additional turbulent viscosity models Were stmulated with values of the mul- 
uply factor MF in equation (3.10) set to 1000, 20000, and 200,000. In addition, a sol- 
ution for the original Hush House model (k- € turbulence model) with the jet engine 


exhaust temperature set 700 ° K was simulated. 
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V. USING THE PHOENICS CODE 


The efficient and successful utihzation of PHOENICS requires a thorough under- 
standing of the PHOENICS Code, the user manual, and the various computing facilities 
required to operate the PHOENICS Code. It 1s extremely difficult to obtain a thorough 
understanding of PHOENICS in a six month period. User manuals in general can not 
provide every bit of information required to fully operate and understand a system. This 
section 1s provided as an additional resource of information to help ehminate common 
meebiems encountered using PFIOENICS. It is intended that this information will re- 


duce the manv setbacks that the beginning users of PHOENICS are sure to stumble on. 


A. ALLOCATING MVS DATA SETS 

The PHOENICS code 1s run as a batch job on the IBM 3033 main frame computer. 
iivesoatch s\stein at the NPS Computing Center 1s MVS (Multiple Virtual System). [he 
MVS system introduces programs to the computer, Initiates their execution, and sched- 
ules all resources and services required. MVS data sets must be allocated prior to run- 
ning the Q! file. For cartesian coordinate system grid, the following data sets are 


required: 
eas o° *  .PHOENTCS.DFO9 
Pmss.5° .PHOENICS.DF10 
Pes 5  .F HOENICS.LOAD 
BSS. FOE NICS SRC 


fete ~~~ should be replaced with the users NPS computer ID number. IVS data 
sets are allocated by use of the CMS command MVSHELP. Select option 2 for allo- 


cation of new data sets. The computer will prompt for the required mformation. If the 


1 


body fitted coordinate (BFC) grid system 1s used instead of the cartesian coordinate grid 


svstem, then two additional data sets are required: 
1. WSS. S244 PL © ue se) alee 
2s NESS: SP Oe SSD 


Table 1 lists the data set parameters for the s:x required files. 


Table 1. DATA SET PARAMETERS 


DATA REGPSI LRECL Bie DSONG 
Set 





Another option of allocating \I1VS data sets and by far the easiest method 1s to re- 
name a previous PHOENICS users data sets using the CMS command MV SITEMPsgr= 
tion 4 the rename option. All MVS data sets are accessible to anvone via the Game 
commands MVSHELP and GETMVS. This procedure will elimmate the problem of 
allocating space for empty data sets for PIHOENICS.DFI2 and PHOENTCSE@eas 
which results in error messages when running the QI file for the first time. 

A few words of caution are needed at this point concerning MVS data sets. If the 
DI-09 file 1s transferred from another user, the restart option in group 24 of the QI file 
should not be used (1.ec..RESTRT should be deactivated) the first time the OY i1leasanam 
unless the data stored on the DFO9 file 1s to be used as input data. After the first 
invoke the restart option. IDDFO9 is the file that the grid geometry and PHOE Nis 
generated output is stored on. Each time the problem ts run the output from the previ- 


ous run 1s read from the DFO09 file as the input data for the present run. The output for 


iow) 
tu 


any run Writes over the previous output on the DF09 file. file. This procedure allows 
sequential runs to be performed with the latest output being submitted as the present 
input data for a restart. In order to save good output files and prevent them from being 
Written over, the DFO9 file must be copied onto a different backup MVS data set. The 
program COPYTAPE hsted in Appendix B provides the means of copving any MVS 
data set to another MVS data set. The new MVS Data File must be allocated prior to 
copying and have the same MVS parameters as the one that is being copied. The name 


of the new file 1s arbitrary but the recommended name form 1s: 
eee. PITOENICS.BACKUP 


Oeitemee ~~ IS the users computer 1D number and BACKUP 1s the new file name (or 
euemmocier ile name). Ihe function of COPYTAPE 1s to copy the SYSUT1 data set to 
mimes 12 data set. COPYTAPE can be used to copy the PHOENICS output file 
Smmeetie backup file or copy the backup file onto the DF09 file which would be used as 
the input data for the next run of the QI] file. It 1s strongly recommended to copy the 
DF09 file onto a backup file on regular intervals such as on a daily or bi-weekly basis. 
It is possible that the DFO9 file could contain bad data as a result of a change in the 
Q1 file or the submission of a different QI file having the same computer ID numbers. 
If this dilemma occurs and the backup files are utilized. onlv a loss of several davs of 


computing would be lost, vice weeks or months if backup files were not uulized. 


B. IMPLEMENTING GROUND 
GROUND is the program for interaction with the EARTH solution program. The 


implementation of GROUND requires five steps: 


1. Insert relevent equations into the file GROUND FORTRAN 
Submit COPYGRD 

Submit COMPGRD 

4. Compact LOAD and SRC MVS data sets 


rt 


Cao 


a3 


5, “Submit iN weer 


Group 9 of the GROUND FORTRAN file 1s shown in Appendix C. Section 5 contains 
the coding required to evaluate ENUT (equation (3.10)) used in the Hush House tur- 
bulence study. The programs COPYGRD, COMPGRD, AND LINKEAR are inel@gied 
in Appendix D for reference. Compacting the LOAD and SRC files requires the se- 
quential submission of three programs for each file. For the LOAD file the following 
prograims listed in Appendix E are required: 

PHOENSOL POADUTE 

PLIOENNOZ EOADEIE 

PFOENNO0s LOADEIB 


For the SRC file the following programs hsted in Appendix E are required: 
PHOBS NOT SOURCERE 
EEOE N02 SOURCE 
EFOE NNO SOONG EEE 
When the LOAD data set has been copied from a previous user it will contain the 
equations placed in the GROUND FORTRAN file by the previous user liga 
equations in GROUND FORTRAN are not required, ensure that the coijimieme 


UCSEGRD= F 1s inserted in group 19 of the Q] file. 


C. DATA TRANSFER FROM MVS TO VAX 

Although the PHOENICS Code 1s available on the NPS computer center ied 
Mainframe, the graphics program PHOTON is only available on the VAX system. 
Therefore in order to have data transfer of the DF09 file from MVS to the VAX 2000 
svstem a data transfer operation is required. Presently there are two methods available 
for data transfer at NPS. The first method involves copying the DFO9 data set toa CMS 
file and then transferring the CMS file via a modem from the mainframe to the VAN 
system. The computer hardware and software for the procedure has not been fully de- 


bugged. More often than not, the transfer of data is interrupted and data transfer 1s 
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seldom successful. Until the svstem has been debugged, it 1s strongly recommended to 
avoid this procedure as itis simply inefficient and undependable. 

The second method available is to transfer the DFO9 MVS data set to a 9 track tape 
and then transfer from the tape to a VAX cartridge. Prior to transferring the DIU9 data 
file to the 9 track tape, the unformatted DI09 file must be formatted so that the file can 
be read by PIIOTON. The program FORMQO9 in Appendix F transforms the original 
unformatted DIF09 file into a new formatted DFO9C file. The new DI 09C file must be 
allocated using MVSHELP prior to submitung FORMO09. ‘The proper setting of the 
parameters for DFO9 are crucial. The data set organization, record format, logical re- 
cord length. and blocksize must be set as listed below. The number of primary and 
secondary storage tracks allocated to the DFO9C file may vary depending on the size of 
(ies O9 file. 

Data set organization = Sequential (PS) 
Record Format = lixed Block (FB) 
Logical record length = 80 

Blocksize = §00 

Primary Allocation = 250 (TRK) 
Secondary = 16(TRK) 

A nine track tape may be checked out from the NPS computer center. A non la- 
beled tape with a density of either 1690 or 6250 bits per inch (bpi) 1s required. If several 
files are to be transfered to tape, the 6250 density tape is recommended duc to its larger 
capacity. The program SEQS5. a uulity program available through the MVSHELP 
command , copies the formatted DI09C file to the 9 track tape. The tape to VAX car- 
tridge transfer equipment is located in the lab of the NPS Meteorology Department 
(Donna Burych is the point of contact for computer access authorization). The follow- 


ing VAX commands are required for data transfer from tape to VAX cartridge: 


CW 
CAN 


Mount, For Record = 80 Blcoksize = 800 MSAO: 
Mount Record = 80 Blocksize = S00 MUAO: 
Copy AIsSAO Sil 3 Ol len 

Dismount MSAO: 

Dismount MUAO: 


During the data transfer ignore the “incompatible attributes” message. It has no effect 
on the transfer operation. When using the cartesian coordinate grid svstem ( which does 
not reqiure the use of MVS DF12 data set ), the DF09 file on the VAX cartridge does 
not require anv editing to be compatible with the Phoenics version 1.4 on the VAX 

In running the PHHOENICS Code it is inevitable that computing errors will occur. 
Listed below are some of the common error codes encountered along with the corrective 
action required: 

Abend code 322 - allocated CPU time has been exceeded. The number of sweeps 
per run must be reduced. In the Hush House model a maximum of 90 sweeps for a 15 
nunute class G job and 360 sweeps for a 60 minute class J job were possible. The time 
per sweep is dependent on the number of equations being solved and the number of it- 
erations of PER tor cache: aniaiele: 

Abend Code SE37 - The LOAD and SRC files require compacting. The compacting 


procedure 1s discussed in section B of this chapter. 
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VI. COMPUTER MODEL RESULTS 


The parametric study of the aerothermodynamics of the NAS Jacksonville, Florida 
Hush House consisted of two independent studies. The first study investigated the ef- 
fects of a 10% reduction of the augmenter tube length on svstem behavior. The second 
study investigated the effects of turbulence model on system behavior. For ease in dif- 
ferentiating the different Ilush House models, the following terms are defined and will 


be referred 1n the discussion of results: 


Model } - Original Hush Flouse Model utilizing the k - ¢ turbulence model (TJET = 
1944 ° K) 


Model 2 - Hush House Model with 10% reduced augmenter tube length utilizing the 
k - ¢ turbulence model (TJET = 1944 ° k) 


Model 3 - Hush House Model utilizing equation (3.10) for turbulent viscosity (MF 
= 200000, TIET = 700 k) 


Nlodel 4 - Hush House Model utilizing equation (3.10) for turbulent viscosity (MF 
= 20000. TJET = 700 k) 


Model 5 - Hush House Model utilizing equation (3.10) for turbulent viscosity (MI 
e000. TJET - 700 k) 


Model 6 - Hush House Model utilizing equation (3.10) for turbulent viscosity (MIF 
= 200, TJET = 700 k) 


Model 7 - Original Hush House Model utilizing k - ¢ turbulence model (TJET - 700 
K) 


The tabulated numerical results generated by PHOENICS for each of the 7 models 
are comprised of nearly 10000 lines of data. For ease of comparing results, output data 
is presented in graphical form for pressure, temperature, and velocity profiles and con- 


EOUTS. 


A. AUGMENTER TUBE LENGTH STUDY 
This study compares the results of model | (the original model) and model 2 (the 


Original model with a shortened augmenter tube). The pressure. temperature, and ve- 
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locity profiles from I1Z= 1-40 at IN=1 and 1Y=4 for model 1 and model 2 are shown 
In Figure 12 on page 43, Figure 13 on page 44 and Figure 14 on page 45. The pressure 
profile shows identical values of pressure for both models from Z=0 to Z= 27.4 meters. 
The pressure in model 2 is approximately 40 pascals higher than in model | from 
Z= 27.4 to Z=45.6 meters. The latter distance represents the 10% shorter augmenter 
tube length. The remaining seetion of the model 2 profile is identical to the model 1 
profile except that it 1s shifted to the left a distance of 2.1 meters, which ts equal to the 
reduction in augmenter tube length. The distance axis is normalized by using the cell 
numbers in the Z direction vice actual distance in the Z direction as Sivemammem 
Figure 15 on page 46. The normalized profile shows more clearly that there is little if 
anv change in pressure magnitudes between the two models. Similarly for the normal- 
ized temperature and velocity profiles shown in Figure 16 on page 47 and Figure 17 on 
page 48 there 1s no significant changes in temperature and Velocity magnitudes between 
the two models. Velocity vector plots, temperature contours, and pressure contours of 
model | and model 2 are shown in Figure 18 on page 49 through Figure 29 on page 
60. Figure 18 on page 49 and Figure 19 on page 50 show the velocity vector plots in the 
test bay at engine level (] Y =4) are identical for model 1 and model 2. Figure 20 om page 
51 and Figure 21 on page 52 show the vector velocity plots at the exhaust Of the 
augmenter tube are identical for model | and model 2. Figure 22 on page 53 and 
Figure 23 on page 54 show the temperature contours in the entrance of the augmenter 
tube are identical for model 1 and model 2. Figure 24 on page 55 and Figure 25 on page 
56 show the temperature contours at the exhaust of the augmenter tube are identical for 
model | and model 2. Figure 26 on page 57 and Figure 27 on page 58 show the pressure 
contours at the entrance of the augmenter tube are identical for model 1 and model 2. 
Figure 28 on page 59 and Figure 29 on page 60 show the pressure contours at the ex- 


haust of the augmenter tube are tdentical for model 1 and model 2. These figures pro- 


Fe 
oe 


i) 


vide verification that the shorter augmenter tube has little effect. on the 
aerothermodynamucs of the flow field. In addition to comparing profiles and contours 
of model | and model 2, the augmentation ratios of each model are compared. The 
augmentation ratio is defined as follows: 

Wee RG 


HIS: : 
Mirve 


where: 


M,. = inlet mass flow rate 


Nigxg = engine mass flow rate 


The inlet mass flow rate is the sum of the mass flow rates through the five air inlet 


Pen) 


Sok 





doors in the test bay. The engine mass flow rate of 38.49 is the same in both 
models. The augmentation ratio 1s a measure of the efficiency of the Hush Ilouse in 
regards to the magnitude of the mass flow rate of secondarv cooling air that can be 
drawn into the Hush House bv the “entrainment effect” of the jet engine exhaust. The 
higher the augmentation ratio the higher the efficiency of the Hush House design. The 
augmentation ratio for model | is 10.66 and for model 2 is 10.87. This amounts to an 
increase of 2°o of the augmentation ratio for model 2 compared to model 1. 

The 10°o shorter augmenter tube is slightly more efficient than the original model. 
Although the pressure. temperature, and velocity of the exhaust gas as it enters the at- 
mosphere is not effected by the shorter tube, the noise level (not analvsed in this thesis) 
is expected to be higher. Recall that the augmenter tube is lined with acoustic absorp- 


tion material. Reducing the augmenter tube reduces the surface area of the absorption 


material thereby reducing the amount of acoustic energy absorbed. 
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B. TURBULENCE STUDY 

The effects of turbulence model on system aerothermal behavior are presented in 
pressure, temperature and velocity profiles and contours. The magnitudes of the turbu- 
lent viscosity, resulting from equation (3.10) , in the arbitrarily selected cell defined bv 
IX=1, 1Y¥=4, and I[Z=20 for model 3, model 4, model 5, model 6, and model 7 are as 


follows: 


model 3. MF= 200,000 726.4 SS 


model 4 MF=20,000 v,=1.1 as 
nc 
pat 
middel 6 L200 enEOen = 
AS 


model 7 k- e model yest ck 


model 5 MF = 1000 Pe etOLONG 


It should be stressed that these values of v,; are representative values of turbulent 
viscosity in one cell onlv. Values of v, are different in each of the cells of the moder 
The results from model 5 and model 6 are essentially the same and are therefore pre- 
sented as a single model referred to as model 5&6. Model 7 (k- ¢ model) 1s included for 
the purpose of providing a relative reference pomt to gage the effects of the temperature 
dependent turbulence models. The pressure profiles for model 3. model 4. model S&6, 
and model 7 are shown in Figure 35 on page 64. In the test bav area prior to the engine 
inlet (Z distance = 0 - 17.6 meters) the pressure profiles for each model are approxti- 
mately constant. The noticable trend 1s that as the turbulent viscosity 1s increased from 
v7=Q0.01 to v;=6.4 the constant pressure profile for each model decreases bv roughly 20 
pascal. The trend continues for the large pressure depression just prior to the engine inlet 
(Z distance = 18.6 meters). After the engine exhaust (Z distance = 22.3 meters). little 
correlation between the pressure profiles and the magnitude of turbulent viscosity is 
observable. For low magnitudes of turbulent viscosity (v;+0.06 and v;=1.1 ) the pres- 


sure profile through the augmenter tube shows little variation. For higher magnitudes 
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of turbulent viscosity (v;=4.4 and v;+6.4 ) the pressure profile through the augmenter 
tube varies from a nagative pressure to a positive pressure. 

The temperature profiles for model 3, model 4, model 5&6. and modei 7 are shown 
in Figure 30 on page 61. The temperature profiles of the four models are identical in 
the test bay (295 ° Kk), through the engine (250 ° kK) to the engine outlet (700 ° kK). It 
is Clearly shown that as the magnitude of turbulent viscosity increases. the temperature 
decrease through the augmenter tube increases as a result of the increased mixing effect. 

The velocity profiles for model 3, model 4, and model 5&6. and model 7 are shown 
in Figure 34 on page 65. The velocity profiles are similar to the temperature profiles in 
that the velocity decrease through the augmenter tube is increases as the magnitude of 
HiewturbuUlent viscosity increases. 

Velocity vector plots in the test bay for model 3, model 4. and model 5&6 are shown 
fevure 35 on page 66, Figure 36 on page 67, and Figure 37 on page 68. The velocity 
vector plot for model 7 1s not included in the figures becuase of its similaritv with the 
meevor plot of model 4. The velocity vectors in model 7 have the same direction but a 
magnitude 1.5 times greater than the velocity vector plots in model 4. 

The vector plot for model 3 shows a smooth flow of air to the engine and augmenter 
tube. The counterflows first observable in Figure 34 ( model 4 ) become more pro- 
nounced in Figure 35 ( model 5&6 ). It appears that in the test bay a higher magnitude 
of turbulent viscosity results in a more distortion free flow path to the engine and 
aucmenter tube. Velocity vector plots at the augmenter tube exhaust for model 3 and 
model 4 are shown in Figure 31 on page 62 and Figure 32 on page 63. The velocity 
vector plots for model 5&6 and model 7 are not shown. Model 5&6 1s identical to model 
4 both in magnitude and direction. The velocity vectors in model 7 are identical in di- 
rection to model 4 but have a magnitude 1|.5 times greater. The gas flow in model 3 exits 


; ra , : a 
to the atmosphere at a velocity of 30 —— and is uniformly distributed across the exhaust 


4] 


opening. The gas flow in model 4 exits at the far end of the exhaust opening at a ve- 
locie Ol a4 — In addition, a counterclockwise vortex is present over the exhaust 
Opening. As in the test bay, a higher magnitude of turbulent viscosity results in a more 
distortion free flow field. 

The temperature contours in the area of the jet engine exhaust for model 3, model 
4, model 5&6, and model 7 are shown in Figure 38 on page 69 through Figure 4] on 
page 72. Temperature contours at the augmenter tube exhaust for model 3, model 4, 
model 5&6, and model 7 are shown in Figure 42 on page 73 through Figure 45 on page 
76. The temperature contours for model 3 (v;=6.4) shows the largest temperature de- 
crease through the augmenter tube ( AT = 382° K ). The temperature decrease through 
the augmenter tube for model 7 (N;=4.4) is A T = 341 ° K. for model 4 (v-=1.1) 1s A 
T = 233° K , and for model 5&6 (v;=0.06) is A T = 173 ° K. The latter eight@iiayiaes 
show that the temperature decrease through the augmenter tube increases as the mag- 
nitude of turbulent viscosity increases. 

The pressure contours in the area of the jet engine exhaust for model 3, model 4, 
model 5&6 and model 7 are shown in Figure 46 on page 77 through Figure 49 on page 
$0. The pressure contours at the augmenter tube exhaust for model 3, model! 4. model 
5&6, and model 7 are shown in Figure 50 on page 81] through Figure 53 on page 84. 
The pressure contours from model 3 and model 7 are similar in that the pressure fields 
change from negative pressures in the augmenter tube entrance to a positive pressure 
at the augmenter tube exhaust. The pressure contours for model 4 and model 5&6 are 
similar in that the pressure fields are negative from the augmenter tube entrance to the 
exhaust stack. It appears that there is a crossover point of the magnitude of the turbu- 


lent viscosity in regards to distribution of the pressure field. For large magnitudes of v; 





NSE : . ae 
(4.4 - 6.4 — ) the pressure changes from negative to positive values. For small magni- 


ANwS . ; 
tudes of v; (0.01 - 1.1] ate the pressure field is always negative. 
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Figure 12. 
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Figure 14. 
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Augmenter Tube Study Pressure Profile (normalized) 


Figure 15. 
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Auginenter Tube Study Temperature Profile (normalized) 


Figure 16. 
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Auginenter Tube Study Velocity Profile (normalized) 


Figure 17. 
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Velocity Vector Plot for Model | (test bay) 


Figure 18. 
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Figure 19. 
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Velocity Vector Plots for model 2 (test bay) 
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Velocity Vector Plot for Model 1 (exhaust) 


Figure 20. 


_ ee a 
| ES A NT 





ale 


cepa rt 8 Ss =a Fe} et 


Pf Soiled AH ot 








eee 


—— eer 


Spe Sa ol PE eae st ee 


et Br 





22 eee ann | A oe ee SP poor a ay es | 


Velocity Vector Plot for Model 2 (exhaust) 
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Temperature Contour for Model | (augmenter tube) 
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Temperature Contour for Model 2 (augmenter tube) 


Figure 23. 
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Temperature Contour for Model I (exhaust) 


Figure 24. 
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Pressure Contour for Model | (augmenter tube) 


Figure 26. 
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Pressure Contour for Model 2 (augmenter tube) 


Figure 27. 
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Pressure Contour for Model 2 (exhaust) 


Figure 29. 
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Turbulence Study Pressure Profiles 


Figure 30. 
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Figure 31. Turbulence Study Temperature Profiles 


62 


O'O0T 





00 


60.0 


S00 


30.0 


DISTANCE(METERS 


(SUA LAW)AONVLSIG 
0°09 0°0S 0°OF o-0e 0°02 o°0T 0°0 


{ 





Goa= 


ereereroevesreeeaneeedg 


OM ie 


eeeseeseeeetoaset®eeee 


O°GL 


year ros ee arate MaeeeatanninTe tht gluta vavcvorevalevotetete\aieisielal sista a sieiale\clelsicialalsleleleiaia\miele ~\a\clwmlala\G.ni@inis\s\eieie)ss\\\es s\e\s 0.4 a\0'*\}\niwielniai® ¥)sis.e:+(a\a\s\c}e s\s\eieiejeieieis.c'c.c ele\nlelssieis sicicisie ce, 0/012)0/irie\~ vine nskcna 


0'Gcl 


« 
whee wee owe cele eee ewe wee ee acess cee eee Pr aetecseseccemseseetesereratesereeseeseegseses 
m 


Oo2T 


0'Sce 


ee eee 8 oe eee ee 
ee ee ——_— e~ ee 


etree eeeeeeeteeeeezeretehoeseereeee eee eateunere s+ Be eave 


O000Z=aN , : | : 
OOQQ002Z=4N : ae Sree way. egies core ere arcroueWeisteresion aire 
CNAOa1 : , : 


0'S2¢ 


YAGNVHO HSNH ATTIANOSMOVS SYN 
V=Al T=XI LV AUAOdd ALIOOTHA 


(S/WDIM ALIOOTHA 


Turbulence Study Velocity Profiles 


a2. 


Figure 


63 


RY Soe RHI 


HHH 


Sy Gl ae seh IS: Aer 





= IT. 





Velocity Vector Plot for Model 3 (test bay) 
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Velocity Vector Plot for Model 4 (test bay) 


Figure 34. 
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Velocity Vector Plot for Model 5 & 6 


Figure 35. 
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Velocity Vector Plots for Model 3 (exhaust) 


Figure 36. 
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Velocity Vector Plot for Model 4 (exhaust) 


Figure 37. 
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Temperature Contours for Model 3 (engine exhaust) 


Figure 38. 


69 


lee ta 


HSMH ADIL edt Ad oe 


Vv |4dCOW 


oS eo ee Se ee ee a ee 


= EA ff 





ems U0 | as (9 


Temperature Contours for Model 4 (engine exhaust) 


Figure 39. 
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Temperature Contours for Model 5 & 6 (engine exhaust) 


Figure 40. 
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Temperature Contours for Model 7 (engine exhaust) 


Figure 41. 
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Temperatuer Contours for Model 3 (exhaust) 


Figure 42. 
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Temperature Contours for Model 4 (exhaust) 


Figure 43, 
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Temperature Contours for Model 5 & 6 (exhaust) 


Figure 44, 
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Temperature Contours for Model 7 (exhaust) 


Figure 45. 
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Pressure Contours for Model 3 (engine exhaust) 


Figure 406. 
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Pressure Contours for Model 4 (engine exhaust) 


Figure 47, 
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Pressure Contours for Model 5 & 6 (engine exhaust) 


Figure 48. 
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Pressure Contours for Model 7 (engine exhaust) 


Figure 49. 
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Pressure Contours for Model 4 (exhaust) 
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Pressure Contours for Model 5&6 (exhaust) 


Figure 52. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 
I. Augmenter Tube Length Study 

The results of the augmenter tube length parametric study shows no significant 
changes of pressure, temperature, or velocity fields when comparing the models with a 
10°o difference in augmenter tube length. System efficiency based on the augmentation 
ratio for the shorter augmenter tube model is 2% higher than in the original model. 
Unfortunately, it is likely that the sound level at the augmenter tube exhaust is inversely 
proportional to the augmenter tube length. That is, decreasing augmenter tube length 
will result in increasing noise levels. Existing jet engine test facilities barely meet, and 
in manv instances exceed, acceptable noise level standards ($5 dba at 200 feet from the 
test cell ). Until a new cost effective acoustic absorption material with a higher coefli- 
cient of absorption is developed or a noise reduction system independent of augmenter 
tube length is developed the design of the Hush Flouse incorporating a shorter 
augmenter tube would be conterproductive. 

2. Turbulence Study 

The selection of a turbulence model as initially expected has a dramatic effect 
on the final flow field solution.As the magnitude of the turbulent viscosity increases, the 
decrease of temperature and velocity through the augmenter tube increase. Halfway 
through the augmenter tube (Z distance = 39.7 meters) model (v;+6.4) has about a 50" o 
reduction in temperature and about a 90% reduction in velocity compared with model 
6 (v;x0.01) which has about a 20° reduction in temperature and about a 35% reduction 
in velocity. The effects of turbulent viscosity on the pressure distribution through the 


augmenter tube are not as clear cut as the effects of turbulent viscosity on temperature 


and velocitv. The magnitude of the turbulent viscositv has a very significant effect on 
the pressure distribution but a simple correlation between the magnitude of the turbulent 
viscosity and the pressure distribution cannot be made. 

It is interesting to note that in general, the results of model 7 (k- € turbulence 
model) fall in a range between the results of model 3 and model 4. The magnitudes of 
the turbulent viscosity ENUT im the arbitrilly selected cell defined by 1X-1, Ty --ama 


IZ = 20 for model 3. model 7. model 4, and model 5 are as follows: 


ENUT for model 3 = 6.4 
ENUT for model 7 = at 
ENUT for model 4 = Lis 
EN itommaeds leone = 


The approximate value of ENUT = 4.4 for model 7 can be arrived at with the multi- 
plving factor MF in equation (3.10) set to the order of 80,000. That 1s the k- e« model 
results could be approxnnated by the temperature dependent turbulence mocelma 
proper setting of MF in equation (3.10). 

This turbulence study was not intended to determine a turbulence model that 
most closely simulated the actual flow field in the Hush House. The study was intended 
to sunply show how different turbulent viscosity models effected the flow field. The 
study shows that to obtain a good numerical solution using PHOENICS, it 1s necessary 
to obtain more experimental data 1n order that the best turbulence model 1s incorporated 


in the PHOENICS simulation. 


B. RECOMMENDATIONS 
The final objective is to develop a valid computer model of the Naval Air Station 
Jacksonville, Florida Hush House. Two major stumbling blocks must be overcome in 


order to attain the final objective. The first and foremost is to obtain more pressure, 


S6 


temperature and velocity data from the augmenter tube of the actual Hlush ITouse. At 
Maeeeml, temperature and velocity data is available for only one location in the 
augmenter tube. The number of locations for collecting the data must be sufficient to 
show accurate pressure, temperature, and velocities profiles through the augmenter tube. 
Once the data is obtained the PHOENICS generated pressure, temperature. and velocity 
profiles should be matched to the actual profiles by manipulating the magnitude of the 
multiply factor MF in the temperature dependent turbulence model. The temperature 
dependent turbulence model is recommended over other turbulence models for future 
use in model simulations of the Hush House due to the relative ease of being able to 


match model results to actual Hush House data. 
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* TAMB: AMBIENT ATR TEMPERATURE "COEGREE oahi iy i) w 

HAMB: AMBIENT ENTHALPY ( JOULES/KG) 7% 

KEINJ: TURBULENT KINETIC ENERGY (JOULES) On sii 

CSUBP: SP HI @ CONST PRE (JOULES, GkG-DEG RELY iNiy) ¥ 
v EPINJ: DISSIPATION RATE OF TURBURENT FE (JCURES OnE we 
¥ KEINA; TURBULENT KINETIC ENERGY (JOULES AMEE * 
w EPINA: DISSIPATION RATE OF TURBULEN® KE (G@U0ULES Jari PE) ees 
* Pl: PRESSURE (NEWTONS/SQUARE METER) # 
¥ Hi: ENTHALPY (JOULES / KG) ve 
* Ul: VEBOCTITY IN X-DIRECMIONSGE eRe, sec) * 
¥ Viv. VELOGITY IN Y-DIRECTIONSCME PERKS, SES) ¥ 
* Wile VELOCITY INSZ4DIREGIION (MNERERS she) " 
¥ TMPi: TEMP OF THE FIRST PHASE VCGDEGR EE Se bin aN) % 
% RHO1: DENSITY OF THE FIRS! PHASE StUkKG /CUBTOSIIE TER * 
“ PRESSO: REFERENCE ATMOSPHERIC “PRESSURES NEY tos = © em ¥ 
we TMPiA: GROUND REF TEMP USED IN GRND2 (DEG KELVIN) 
* TMPIB: 1/CSUBP USED IN GRNDZ EOQU (KG-TPEG TE LY IN yours ¥ 
¥ DRHIDP: DELTA RHO/DELTA PRESSURE W,/ RESPECT iam % 

FACT: CORRECTION OR "FUDGE" FACTOR * 
os ENUL: LAMINAR KINEMATIC VISCOSITRY (MEDERS S0/SE@) w 
E Nae F eaearaan eee VISCOSITY Cee ee 
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GROUP 1: KUN TITLE AND OTHERS PREC EMINAK TES 

TEXTCNAS JACKSONVILLE HUSH CHAMBER) 

REAL(WJET ,MJET , TJET , PJET , RHOJET , HJET , KEJET , GASCON , CSUBP,, TAMB , HAMB) 

REALCKEINJ ,EPINJ,KEINA,EPINA) 


GROUP 2. TRANSIENCE; TIME-STEP SPECIFICATION (DEPAULI=STEADY) 
GROUP 3. X-=DIRECTION GEID SPECICATION 
NX=10 


XFRAC(1)=-1 ; XFRAC(2)=0. 31 
XFRAC(3)=4 ; XFRAC(4)=0. 5475 


&$ 


XFRAC(5)=3 3; XFRAC(6)=2. : 
XFRAC(7)=2 3; XFRAC(8)=1. 


GROUP 4. Y-DIRECTION GRID SPECIFICATION 


NY=12 
Meee i j=-3 <; YFRAC(2)=0. 562 
YFRAC( 3 )=2 | YPRACC4)=0, 314 
YFRAC(5)=3 ; YFRAC(6)=0. 562 
YFRAC(7)=4 ; YFRAC(8)=0. 75 
GROUP 5. Z-DIRECTION GRID SPECIFICATION 
NZ=40 
ZFRAC(1)=-1 3; ZFRAC(2)=1.0 
ZFRAC(3)=1 ; ZFRAC(4)=1. 7 
ZFRAC(5)=1 ; ZFRAC(6)=2. 0 
ZFRAC(7)=1 ; ZFRAC(8)=1. 7 
ZFRAC(9)=2  ; ZFRAC(10)=3.0 
ZFRAC(11)=1 ; ZFRAC(1i2)=1.5 
ZFRAC(13)=1 ; ZFRAC(14)=1. 0 
ZFRAC(15)=1 3; ZFRAC(16)=1. 7 
yepeeeCl7)=2 ; ZFRAC(18)=1. 0 
ZFRAC(19)=1 ; ZFRAC(20)=1. 7 
memAGC(Z1)=2 ; ZFRAC(22)=1.0 
ZBRAC(23)=3 ; ZFRAC(24)=0. 566666 
ZFRAC(25)=2 3; ZFRAC(26)=0. 35 
ZFRAC(27)=7 ; ZFRAC(28)=3. 00 
ZFRAC(29)=1 3; ZFRAC(30)=2. 0 
Peeec(31)=4 ; ZFRAC(32)=0.75 
ZFRAC(33)=4 ; ZFRAC(34)=0.575 
Peee(35)=4 ; ZFRAC(36)=0.5 
ZFRAC(37)=1 3; ZFRAC(38)=3. 0 


GROUP 6. BODY-FITTED COORDINATES OR GRID DISTORTION 


GROUP 7. VARIABLES STORED, SOLVED & NAMED 
ome (CP 1 ,H1,U1,V1,W1) 
SOLUTN(P1,Y,Y,Y,N.N,N) 
eee 2MP1,RHO1,ENUT) 


GROUP 8. TERMS (IN DIFFERENTIAL EQUATIONS) & DEVICES 


miaoUP 9. PROPERTIES OF THE MEDIUM 
meMB=295. 
TJET=700 
PJET=0. 
P= 19/7. 72 
Pepeo0—101325. 
RHO1A=0. 
GASCON=286. 
RHO1B=1. /GASCON 
fme. & 1=(( PRESSO+PJET )**RHOIB ) /TJET 
PET =MIET /RHOJET 
TMP1A=TINY 
~=- TIMPiB IS EQUAL TO 1. DIVIDED BY C-SUB-P; USED IN TMP1=GRND5 
CSUBP=1004. 
HAMB=TAMB*CSUBP 
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TMP1B=1. /CSUBP 
HJET=CSUBP*TJET 
wee THE NEXT 2 STATEMENTS INVOKE COMPRESSIBLE FLOW IN GREXK1 
TMP 1=GRND2 
RHO1=GRND5 
DRH1DP=GRND5 
ENUT=GRND 
KEINJ=0. 5*(0. OOS*WJET)**2 | EPINJ=0.09" REIN? 57 (005-0 
KEINA=0. 5*(. 005*1. 0)**2 ; EPINA=0. 09*KEINA« 1) 5/7 G@G=a) 


GROUP 10. INTER-PHASE-TRANSFERW PROCESSES AND VP R@r angie 


GROUP 11. INITIALIZATION OF VARIABLE OR POROSITY FIELDS 
FIINIT(H1)=HAMB 
FIINIT(TMP1)=TAMB ; FIINIT(DEN1)=1. 0 
FIINIT(KE)=KEINA ; FIINIT(EP)=EPINA 
FIINIT(P1)=0. 0 
FIINIT(W1)=0. 0 
FIINIT(V1)=0. 0 
FIINIT(U1)=0. 0 
PATCH(WINIT, INIVAL,1,5,1,8,15,30,1,1) 
INIT(WINIT,W1,0. 0,0. 2*WJET) 
INIT(WINIT,H1,0. 0,CSUBP*0. 2*TJET) 
INIT(WINIT,KE,0. 0,KEINJ*0. 2) 
INIT(WINIT,EP,0. 0, EPINJ*0. 2) 
wets SEE GROUP 13 FOR POROSITY FIELDS 


GROUP 12. CONVECTION AND DIFFUSION ADJUSTMENTS 


GROUP 13. BOUNDARY CONDITIONS AND SPECIAL SOURCES 
DEFINE SOURCES AND SINKS FOR JET ENGINE? 
CONPOR(0. 0,CELL,1,1,4,5,12,14) 
PATCH( JETIN HIGH) 11-205 e eres) 
COVAL( JETIN,P1,FIXFLU, -MJET) 
COVAL( JETIN,H1,ONLYMS , SAME) 
PATCH( JETOUT, LOW, 121-45 15 steele) 
COVAL( JETOUT ,P1,FIXFLU ,MJET) 
COVAL( JETOUT ,H1,ONLYMS , HJET) 
COVAL( JETOUT ,W1,ONLYMS ,WJET) 
COVAL( JETOUT, KE ,ONLYMS , KEINJ) 
COVAL( JETOUT, EP , ONLYMS , EPINJ) 


“EAST BOUNDARY TO ATMOSPHERE? 
PATCH( OUTSIDE. BAST 10, 10 ite 201 mae 
COVAL(OUTSIDE,P1,0.1,0. 0) 

COVAL( OUTSIDE ,H1,FIXVAL, HAMB) 


"TOP OF HUSH HOUSE AND TUBE INCLUDING TUBE LIP** 

CONPOR( 0. 0,NORTH,1,9,12,12,1,19) 

CONPOR(0. 0,NORTH,1,5,-8,-8,18,27) 
CONPOR( 0. 0,NORTH,1,5,-2,-2,35,35) 
CONPOR( 0. 0,NORTH,1,5,-4,-4,36, 36) 
CONPOR( 0. 0,NORTH,1,5,-5,-5,37,37) 
CONPOR(0. 0,NORTH,1,5,-6,-6,38, 38) 
CONPOR(0. 0,NORTH,1,5,-7,-7,39, 39) 
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mw 2OP OF MODEL TO ATMOSPHERE** 
fewer 1 OPOUT,NORTH,1,10,12,12,20,40,1,1) 
Bev AE CLOPOUT,P1,0.1,0.0) 
COVAL( TOPOUT , H1 ,ONLYMS , HAMB) 
eee TOPP,NORTH,10,10,12,12,1,19,1,1) 
eMac  TOPP,P1,0.1,0.0) 
COVAL(C TOPP ,H1,ONLYMS , HAMB ) 


**FLOOR OF HUSH HOUSE, TUBE AND GROUND** 
oe OR 0. 0,SOUTH,1,10,1,-1,1,40) 


**FORE AND AFT BULKHEADS IN HUSH HOUSE AND TUBE** 
CONPOR(0.0,LOW,1,9,1,12,-1,-1) 
CONPOR( 0. 0,LOW,6,9,1,12,-20,-20) 
CONPOR(0. 0,LOW,1,5,9,12,-20, -20) 
CONPOR(0.0,LOW,1,5,1,2,-35,-35) 
CONPOR(0.0,LOW,1,5,3,4,-36,-36) 
MONPOR(0.0,LOW,1,5,5,5,-37,-37) 
CONPOR(0. 0,LOW,1,5,6,6,-38, -38) 
CONPOR(0.0,LOW,1,5,7,7,°39,739) 
CONPOR( 0. 0,HIGH,1,5,8,8,-39,-39) 


‘©VERY FRONT MODEL STRIP TO ATMOSPHERE" 
PATCH( FRONT, LOW,10,10,1,12,1,1,1,1) 
COVAL( FRONT ,P1,0.1,0.0) 

COVAL( FRONT ,H1,ONLYMS , HAMB) 


**THE VERY END OF OUR MODEL TO ATMOSPHERE** 
Peden ENDEND,HIGH,1,10,1,12,40,40,1,1) 
ee b(ENDEND,P1,0.1,0.0) 
COVALCENDEND,H1,ONLYMS , HAMB ) 


**OUTSIDE EAST WALL IN HUSH HOUSE** 
meee 0.0,EAST,-9,-9,9,12,1,19) 
BeeeORCO.0,EAST,-9,-9,1,8,1,1) 


CONPOR(0. 0,EAST,-9,-9,1,8,3,3) 
CONPOR( 0. 0,EAST,-9,-9,1,8,5,8) 

CONPOR(0. 0,EAST,-9,-9,1,8,10,11) 
CONPOR(0. 0,EAST,-9,-9,1,8,13,14) 
CONPOR( 0. 0,EAST,-9,-9,1,8,18,19) 


OUTSIDE OF TUBE PLUS TUBE LIP IN HUSH HOUSE** 
memeor(0.0,EAST,-5,-5,1,8,18,39) 


GROUP 14. DOWNSTREAM PRESSURE FOR PARAB=. TRUE. 


BEOUP 15. TERMINATION OF SWEEPS 
FSWEEP=1 ; LSWEEP=90 


GROUP 16. TERMINATION OF ITERATIONS 
LITER(P1)=40 
GROUP 17. UNDER-RELAXATION DEVICES 


REAL( DELT, FACT) 
DELT=50. 0/(NZ*0. 2*WJET) 
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FACT=10. 
FACT=1. 00 
FACT=0. 7 
FACT=3. 5 
FACT=5. 
FACT=0. 50 
FACT=0. 1 
RECAXGE Pe oINR EX Cue 
RELAX(U1,FALSDT , FACT*DELT ) 
RELAX(V1,FALSDT , FACT*DELT ) 
RELAX(W1, FALSDT , FACT**DELT ) 
RELAX(C HI, FPALSDT PACT iin) 
RELAXCKE , FALSDT , FACT**DELT) 
RELAX( EP, FALSDT , FACT*DELT ) 


S 


GROUP 18. LIMITS ON VARIABLES OR INCREMENTS TO THEM 
VARMIN(U1)=-300. 
VARMAX(U1)=1000. 
VARMIN(V1)=-300. 
VARMAX(V1)=1000. 
VARMIN(W1)=-100. 
VARMAX(W1)=1500. 

VARMIN( TMP1)=250 
VARMAX( TMP1)=TJET 
VARMIN(H1)=1000. 
VARMAX(H1)=CSUBP*TJET 
VARMIN(P1)=-3. 5E+05 
VARMAX(P1)= 3. 5E+05 
VARMIN(RHO1)=. 05 
VARMAX(RHO1)=S. 
VARMIN(KE)=1. E-10 
VARMAX(KE)=1. E10 
VARMIN(EP)=1. E-08 
VARMIN( VIST)=1. E-20 
VARMAX(EP)=1. 3**(KEJET*#*1. 5) 


GROUP 19. DATA COMMUNICATED BY SATELLITE TO GROUND 
USEGRD=T 


GROUP’ 20. PRELIMINARY PRENM SOUT 
GROUP 21. PRINT-OUT OFeVARIABLES 


GROUP 22. SPOT-VALUE PRINT-OUT 
“MONITOR THE OUTPUT AT DESIGNATED POINTS?" 
IXMON=1 ; IYMON=4 ; IZMON=16 
TSTSWP=1 ; LUPR3=6 
ABSIZ=1.0 ; ORSIZ=1. 0 


GROUP 23. FIELD PRINT-OUT AND PLOT CONTROL 
NUMC LS=10**NXPRIN 
ITABL=3 ; NPLT=1 ; NZPRIN=1 
IPLTL=LSWEEP 
NPRINT=3 


GROUP 24. DUMPS FOR RESTARTS 
SAVE=T 
feo ORECP1,U1,V1,W1) 
RESTRT( ALL) 
SLOP 


// 
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APPENDIX B.) PROGRAM FOR COPYING MVS DATA SETS 
A. COPYTAPE 


//RATNER JOB (2529,9999),'XHUSH' ,CLASS=A 

/ | CXESEC PGM=IEBGENER 

//SYSPRINT DD SYSOUT=A 

//SYSIN DD DUMMY 

//SYSUT1 DD DSN=MSS. 51541. PHOENICS. CONV7 , DISP=SHR 

| FavVsoi2 DD DSN=MSS.51541. PHOENICS. DFO9G ,DISP=( SHR) 


fe 


94 


APPENDIX C. “GROUND? FILE (GROUP 9) 


Section 35 ON > ae oes ste see to ev gues Do Oe (eqation (3,10): 
C* rete irint ete te ele tire ibicibini ee eile rie bit teri brn tt he) he peter a tela bil kta lta rit tari riick ata tk vob titoket deh toni 
© 
G=—=- GROUP 9. Properties of the medium (or media) 


C 
C The sections in this group are arranged sequentially in their 
C order of calling from EARTH. Thus, as can be seen from below, 
C the temperature sections (10 and 11) precede the density 
C secetons (1 and 3); so, density formulae can refer to 
C temperature stores already set. 
- oe TO Cae os ee ne nee 96 ad rec oan a » 900, Se ee ieee ISC 
C rae aeICs! Rice riety eter iieicl eo lsiclanteeinr rise or x bevy Rese Steecs ook Fgeltche tke sels Seis 
900 “CONTINUE 
ee ee 2 SECTION 10 -----eer crc rer errr reer rere ee- 
C Poet LeaGkRND---———=—— phase-1l temperature Index AUX(TEMP1) 
RETURN 
901 CONTINUE 
eee ~~ ~~~ ee SECTION 11 errcecn ccc rer r cer err ere---- 
6 For TMP2. LE. GRND--------- phase-2 temperature Index AUX(TEMP2) 
RETURN 
902 CONTINUE 
ees SbCl 2 ===. = = = 
C Boast LE. GRND--=—------— phase-1 length scale Index AUX(LEN1) 
RETURN 
ae CONTINUE 
=  —  — SHOT ROIS, UG OC SS 6 6 i tS 
C For Ib) We, (Ee SS or phase-2 length scale Index AUX(LEN2) 
RETURN 
91 CONTINUE 
fees ooo moon m ee ---- SECTION 1 <"9-333< renee nnn------=---- 
C For RHO1. LE.GRND--- density for phase 1 Index AUX(DEN1). 
965 CONTINUE 
RETURN 
92 CONTINUE 
Reo oo em eee ee eee e ee SIO UIIIIOIN eS NN a 
G For DRH1DP. LE. GRND--- DCLN(DEN))/DP for phase 1 (DI1DP). 
RETURN 
Poe CONTINUE 
i Mites a ie ee SE LON my ae = <a ae = l= = Sie = + = =~ 
C For RHO2. LE. GRND--- density for phase 2 Index AUX(DEN2). 
RETURN 
94 CONTINUE 
GC * een ncern cc cnner-e---- DEAes tO) Cee = = im a le oe ie i i 
6 For DRH2DP. LE.GRND--- DCLNCDEN))/DP for phase 2 (D2DP). 
RETURN 
95 CONTINUE 
GF eecr rece reece rene -- SSA 9) ae ene ee = 
C For ENUT. LE.GRND--- reference turbulent kinematic viscosity. 


Combe noon Viol) ,AUXC TEMP) 54. 666,5.432E-2 "2. 426E-5,7.93i1E-9, 


DS 


1-1. 104E-12) 

CALL FN25(AUX(VIST),1. OE-6) 
CALL FN27( AUX(VIST) , AUX(DEN1) ) 
CALL FN25(AUX(VIST) , 1000. 0) 
RETURN 


wee rence eee e------- SEGl 1 ON 6 ee 
For ENUL. LE. GRND--- reference laminar kinematic viscosity. 
RETURN 
97 CONTINUE 
I Te Oe ea ee 
For PRNDTLC ).LE.GRND--- laminar PRANDITL nos. > er diffusive 
RETURN 
98 CONTINUE 
9 SIG@INUO Mg) Season seo See See SS See Ses So sSe 
For PHINT( ). LE.GRND--- interface value of first phase(FII1). 
RETURN 
99 CONTINUE 
HS SSS SSS SS SSeS SS SEGDION | 90 =o a= eo i 
For PHINT( ). LE.GRND--- interface value of second phase(FII2) 
RETURN 


eis 


Crete dedede dere ve te dete dete de Te Te Fe Fe Fe Fe Ve Fee Be Ve Te ETE TE TE TE TE TE TE TE TS VOTE Te TE FE VOTE TE TENE ETE TE TE TEE TE TE TS PEGE PEGE GE BEET TE TE 
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APPENDIX D. PROGRAMS REQUIRED FOR IMPLEMENTING 


GROUND 


A. COPYGRD 
//COPYGRD JOB (2529,9999),'COPYTMVS' ,CLASS=A 
// EXEC PGM=IEBGENER 
//SYSPRINT DD SYSOUT=A 
//SYSIN DD DUMMY 
//SYSUT2 DD UNIT=SYSDA, 
DSN=MSS. $1541. PHOENICS. SRC( GROUND) , DISP=SHR 
PSYSUT1 DD * 
PROGRAM MAIN 
C THIS IS THE MAIN PROGRAM OF EARTH 
MEBEARENAME GROUND. PIN-=--=------=-=-24-.---42-----_-- 16 July 1986 


~~ 
~~ 


mo COPYRIGHT 1984, LAST REVISION 1986. 

CONCENTRATION HEAT AND MOMENTUM LTD. ALL RIGHTS RESERVED. 
This subroutine and the remainder of the PHOENICS code are 
proprietary software owned by Concentration Heat and Momentum 
Limited, 40 High Street, Wimbledon, London SW19 5AU, England. 


PROGRAM MAIN 


1 The following two COMMON's, which appear identically in the 
satellite MAIN program, allow up to 25 dependent variables to 
be solved for (or their storage spaces to be occupied by 
other variables, such as density). If a larger number is 
required, the 25's should be replaced, in the next 8 lines, 
by the required larger number; and the 100 in COMMON/FO1/ 
should be replaced by 4 times the required number. Numbers 
less than 25 are not permitted. 


CC) -GasC Gata a Ce Ca Ga Ca a Cea 


COMMON/LGE1/L1(25) /LGE2/L2( 25) /LGE3/L3(25)/LGE4/L4(25) 
1/LDB1/L5(25)/IDA1/11(25) /IDA2/12( 25)/IDA3/13(25)/IDA4/14( 25) 
1/IDA5/15(25) /IDA6/16(25)/GI1/17(25)/GI2/18(25)/HDA1/IH1(25) 
1/GH1/1H2( 25) /RDA1/R1( 25) /RDA2/R2( 25) /RDA3/R3( 25) /RDA4/R4( 25) 
1/RDA5/R5( 25) /RDA6/R6( 25) /RDA7/R7(25)/RDA8/R8(25)/RDA9/R9(25) 
1/RDA10/R10( 25) /RDA11/R11(25) 
1/GR1/R12( 25) /GR2/R13( 25) /GR3/R14(25)/GR4/R15(25) 
1/IPIP1/1P1( 25) /HPIP2/1HP2( 25) /RPIP1/RVAL( 25) /LPIP1/LVAL( 25) 
1/IFPL/IPLO( 25) /RFPL1/ORPRIN( 25) /RFPL2/ORMAX(25) 
1/RFPL3/ORMIN( 25) /RFPL4/CELAV(25 ) 

ROCIO AM Wilen2 13, bo iS ePBChI nL VAL 

CHARACTER*4 IH1,IH2,IHP2,NSDA 


C 
COMMON /F01/19( 100) 
COMMON/DISC/DBGFIL 
EXTERNAL WAYOUT 

C 


on 


Grz Set dimensions of data-for-GROUND arrays here. WARNING: the 
C corresponding arrays in the MAIN program of the satellite 
g (see SATLIT) must have the same dimensions. 
COMMON/LGRND/LG(20)/IGRND/IG(20)/RGRND/RG( 100) /CGRND/CG( 10) 
LOGICAL LG 
CHARACTER*4 CG 
C 
Ga Set dimensions of data-for-GREX1 arrays here. WARNING: the 
C corresponding arrays in the MAIN program of the satellite 
e (see SATLIT) must have the same dimensions. 
COMMON/ LSG/LSGD( 20) /ISG/ISGD( 20) /RSG/RSGD( 100) /CSG/CSGDC 10) 
LOGICAL LSGD 
CHARACTER*4 CSGD 
C 
C 4 Set dimension of patch-name array here. WARNING: the array 
C NAMPAT in the MAIN program of the satellite must have the 
C SALE DIMENSION: 
COMMON/NPAT/NAMPAT( 500) 
CHARACTER*8 NAMPAT 
C 
C Declare local CHARACTER variables. 
CHARACTER NDUM4*4 ,NDUM6*6 , NDUM15*15 
C 
Ge The numbers in the next two statements (which must be ident- 
C ical) indicate how much computer memory is to be set aside 
C for storing the main and auxiliary variables. The user may 
© alter them if he wishes, to accord with the number of 
C grid nodes and dependent variables he is concerned with. 


COMMON F(500000) 
NFDIN=500000 


Cro Logical-unit numbers and file names, not to be changed. 
DBGFIL=. FALSE. 
CALL DSGEAR(14,LUPR3, . 1oeNDULMG ieewiG) 
CALL DSCEAR(6,LUDUM,' ',4,NDUM4,9,33) 
CALL DSCEAR(-10,LUSDA,' ',4,NSDA,0,0) 
CALL DSGEAR(-14,LUPRI, 745 NDNI5..0. 0) 
CALL DSGEAR(21,LUDST,° ~ ,4,NDUM4 ose) 


User may here change message transmitted to logical unit 

LUPR3 

CALL WRIT40('GROUND STATION IS GROUND. FTN 11 JULY 86 '‘) 

CALL MAINICNFDIM, LUPR1, LUPR3 , LUSDA,NSDA) 

CALL WAYOUT(0) 

S10 

END 

Ce Pereve te ee Tete TT ETE TS TEE TE TEE TE TE TS TS TE TEE TE TS Te FETE TE VE FETED TETE TE TENE PE TENE TC TS TE TEES ECC TE TE PETE TET TEE TEE TE TEE TE 

SUBROUTINE GROSTA 

€ (C) COPYRIGHT 19846, LAST REVism@en iced 

C CONCENTRATION HEAT AND MOMENTUM LTD. ALL RIGHTS RESERVED. 

INCLUDE (SATEAR) 

INCLUDE (GRDLOC) 

INCLUDE (GRDEAR) 

This subroutine directs control to the GROUNDs selected by 

the satellite settings of USEGRX, NAMGRD & USEGRD. 
Subroutine GREX1 contains much standard material, eg. 


Co Gaca 


<2 0) 
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Ce) Go) Sep Hee, Gp a Gp zs >} Dp ie es Sp Cl) CG OG) amaaa Co 020) 


options for fluid properties, several turbulence models, 
wall functions, etc. 


IFC USEGRX) CALL GREX1 


ESTER is for electrolytic-smelter modelling of the Hall-cell 
and Soderberg types used in the reduction of aluminium. 


IF(NAMGRD. EQ. 'ESTR') CALL ESTRGR 


SCRS contains the simple-chemical-reaction-model of 
combustion, the theoretical basis of which is found in the 
book "Combustion & Mass Transfer" by D B Spalding (1979) 
This ground also contains geometrical features of a 
simplified can combustor. 


IF(NAMGRD. EQ. 'SCRS') CALL SCRSGR 
A more advanced model of a combustor is given in COMBGR. 
IF(NAMGRD. EQ. 'COMB') CALL COMBGR 


WJETGR shows how to represent non-isotropic effects in the 
turbulence of a wall jet. 


IF(NAMGRD. EQ. 'WJET') CALL WJETGR 


TRACGR contains software for tracking fluid interfaces by 
means of a set of imaginary particles which follow the motion 


IF( NAMGRD. EQ. 'TRAC') CALL TRACGR 


PARTGR is used to solve for the motion of particles slipping 
relative to the host fluid. A spectrum of particle sizes can 
can be represented. Each particle is characterized by a size, 
an interphase friction coefficient, an evaporation rate & a 
temperature. 


IF(NAMGRD. EQ. 'PART') CALL PARTGR 


RADIGR provides the coding sequences required to activate 
the so-called six-flux radiation model. 


IF(NAMGRD. EQ. 'RADI') CALL RADIGR 


GAUSGR provides the Gauss-Seidel solver as an alternative 
to the whole-field linear equation solver provided in EARTH. 


IF(NAMGRD. EQ. 'GAUS') CALL GAUSGR 
NOZLGR provides initial conditions & special print out for 
a convergent-divergent nozzle case for which body-fitted 


coordinates are used. 


IF(NAMGRD. EQ. 'NOZL') CALL NOZLGR 


9? 


ye AEROGR provides inlet boundary conditions & initial conditions 
C for a one-half C grid for dn aeroroii, 
C 
IF(NAMGRD. EQ. 'AERO') CALL AEROGR 
C 
Cas POLRGR specifies uniform flow boundary conditions into 
C a polar domain of 360 degree extent. 
C 
IF(NAMGRD. EQ. 'POLR' ) CALL POLRGR 
C 
ce BTSTGR contains the sequenses used in conjunction with 
C the BFC test battery. 
C 
IF(NAMGRD. EQ. ‘BTST')} CALL BTSTGR 
C 
C2. TESTGR contains test battery sequences used in conjunction 
C with the test-battery SAILII subroutine, [FESTS!. 
C 
IF(NAMGRD. EQ. 'TEST') CALL TESTGR 
C 
Cae SPECGR is a generic special’ GROUND the name of which can 
C be used by anyone for their own purposes. 
C 
IF(NAMGRD. EQ. 'SPEC') CALL SPECGR 
C 
Ca. The model ground is for the insertion of new user sequences. 
C 
IF(USEGRD) CALL GROUND 
C 
Ca The data echo is now called at the preliminary print stage. 
C 
IFC IGR. NE. 20) RETURN 
TLE CNOT ECHO mCOs1 Os 
CALL: DATPRNCY 3 Y5¥ 505, 5 Yop We ete eye 
1 XV g \eoraer At Weenie) 
RETURN 
20 CALL DATPRNCY,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N,N) 
RETURN 
END 
Cee se sete tess Tee Bete STE TE TE TE TS ETE TSS ESTEE TOTS ISTE TE TE TESTE TE TE ESTE TE TE TS TS TE TE TS TE TE TESTE TS TE TET AEE TE TESTS TEBE TET ESTE 


SUBROUTINE GROUND 
C (C) COPYRIGHT 1984, LAST REVISTON=193e. 
C CONCENTRATION HEAT AND MOMENTUM LTD. ALL RIGHTS RESERVED. 
INCLUDE (SATEAR) 
INCLUDE (GRDLOC) 
INCLUDE (GRDEAR) 
INTEGER HIGH,OLD,AUX,SOUTH,NORTH,EAST,WEST 
LOGICAL STORE ,SOLVESPRING 
CXXXXXXXXRAXXAXXAXXAXXXXXXXXXXXAXXXAXXKXX USER SECTION STARTS: 


1 Set dimensions of data-for-GROUND arrays here. WARNING: the 
corresponding arrays in the MAIN program of the satellite 
and EARTH must have the same dimensions. 
COMMON/LGRND/LG( 20) /IGRND/IG(20)/RGRND/RG( 100) /CGRND/CG( 10) 
LOGICAL LG 


C)-C).C) Co 


100 





CHARACTER*4 CG 


C 

C2 User dimensions own arrays here, for example: 

C Pie ENelON UUHC1O,10),UUC(10,10).,UUXC 10,10), UUZC 10) 
PARAMETER( MY=20 , MX=1) 
DIMENSION GVICMY ,MX) ,GW1CMY,MX) ,GTMP1( MY ,MX) ,GRHO1(MY,MX) 
DIMENSION GEICMY,MX) ,GENULCMY ,MX) ,GMACH( MY ,MX) ,GHCCMY,MX) 
Diba tON GSTCMY,MxX),GCFCMY ,MxX) ,GOCMY ,MX) ,GTAUC MY , MX) 
COMMON CP,GAMA,RAIR, PRL, RFAC 


3. User places his data statements here, for example: 
DATA NXDIM,NYDIM/10,10/ 


C2 Ci Gort) Co 


4 Index functions for GROUND-EARTH variable references. 
LOW(C I )=NPHI+1 
HIGH( I )=2**NPHI+I 
OLD( 1)=3*NPHI+! 

INC 1 )=4**NPHI+I 

STORE( I)=MODCISLNC I) ,2). EQ. 0 
SOLVEC L)=MODCISLNC 1) ,3): EQ. 0 
Pee PGi y= 00( 1PEN( 1), 2), E00 
SOUTH( I)=-(KFCI)-1) 
NORTH( I)=-C KFC 1)+4+1) 

EAST( I)=-(KFC1I)+NY) 

WEST(C I)=-(KFCI)-NY) 


5 Insert own coding below as desired, guided by GREX1 examples. 
Note that the satellite-to-GREX1 special data in the labelled 
COMMONs /RSG/, /1ISG/, /LSG/ and /CSG/ can be included and 
used below but the user must check GREX1] for any conflicting 
uses. The same comment applies to the EARTH-spare working 
Pemeacme oot ee bior? ,. 7. HAoP MO, If the call to GKEX] has been 
deactivated then they can all be used without reservation. 


C2 GG) - 6) Gl te) 


IXL=IABS( IXL) 
fPOoIGk EO. 13) GO TO 13 
ek EO eto) 160 TO 19 
MMC nse G05 6.7.6.9 10.1) 19. 13 114 15,196.17, 18,19,20,21, 
lee, 923 9 2a) IGR 
Cvevevededede deve ceve ds vere we 7 Sehr iris tartare oo i Gr iahirtis lr orGicnracbir tiara berber brrarhiie oils cuir hark far il 
C 
C--- GROUP 1. Run title and other preliminaries 
C 
mecen lo 1001,1002),1SC 
1001 CONTINUE 


RETURN 
1002 CONTINUE 
RETURN 
Ciedededededs teds de de deve ede fev ve de ve ve Fee FeFe TE Te TENS Ve TENE TE TE TENE TET HEED TELE TEE TE TE TE FETE TE TE TENS ETE ETE EVE VEN DETECTS 
C 
C--- GROUP 2. Transience; time-step specification 
C 
2 CONTINUE 
eee 
Casdedvde deve decks deve de decode de deve ds ve deve de te Fev Fe Te Fe Ve Ne Fe FE FE TENE TE TET TE TERE TE TE TE TE TE VE FE FE TENE NE TE TE TEE TELE TE TE TE EE HE 
C 
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C--- GROUP 3. X-direction grid specification 


G 
3 CONTINUE 
eae 
Crvdedededevedsvevede deve vedeie He %e% Kesleseskvovvevevesedevesevdesecdcsesevcslesesleokstevievesivestcocsesese seve okeseslesieslcsicsiesieskotesksesc 
C 
C--- GROUP 4. Y-direction sridespecisicaeien 
C 
4 CONTINUE 
RETURN 
CO Bere tess tess ds Tess Tove Ve Te Te Te Te Te Te TEENS FETE TE Fe TE Fe To ETE TE TEE TE Te TE TE Me VEE TENE TE He ETE VEE TS TENE TE GETS Ve FETE Fe Ne FFE Tee TE 
C 
C--- GROUPS. Z2-dareceion orid) specwricaeaon 
C 
5 CONTINUE 
RETURN 
Crevidevedescasdededs devededs veut as ve ved ae ve ve ve ve deve de He deve deve ve weve He seve Jove ve deve Heys vewe ve ve seve TeNe Fe Tete Te dese vege ae se ve 
c 


C--- GROUP 6. Body-fitted coordinates or grid distortion 
C 
6 CONTINUE 
an 
C wevevevedeveds dodevs seve ve de ve de de de vedede deve ds Fe ve ve ve ve Tee ve ve seve Te ve se de ve seve Tose Ue Fe Vee Fe Fe Vege Fe Fe Fe eee Fe Fe eve ye Fs He 
C 
C--- GROUP 7. Variables stored, solved & named 
7 CO xEINUE 
RE PFURN 
Crete sede Tee Tete Se Te Te TET TEE TEETER EEE EEE BBE EEE EBA ERE EEE REE BER BEBE DERE EERE 
C 
C--- GROUP 8. Terms (in differential equations) & devices 
G 
8 GO TO (81,82,83,84,85,86,87,88,89,810,811,812,813, 814,815) 
lise 
a CONTINUE 


Co ele im eee ei SECTION 1 rece cr enc c csc ccscrceeren--- 
C Slo U1AD. LE.GRND--- phase 1 additional velocity (VELAD). 
RETURN 
= CONTINUE 
Cae SECTION 2 <=<<<--------- -— aa 
C ee U2AD. LE. GRND--- phase 2 additional velocity (VELAD). 
RETURN 
83 CONTINGE 
Ge Sa oo ae SC UO) a 
C For V1iAD. LE. GRND--- phase 1 additional velocity (VELAD). 
RETURN 
84 CONTINUE 
OMS a SS on oeaae SECTION 4 crccerrrrercereee-------=2- 
S For V2AD. LE. GRND--- phase 2 additional velocity (VELAD). 
RETURN 
85 CONTINUE 
Gee Soa poe assSoSoscsosess SC LON i ee 
C For WIAD. LE.GRND--- phase 1 additional velocity (VELAD). 
RETURN 
86 CONTINUE 
OM OS SECTION 6 === - ee 


C For W2AD. LE. GRND--- phase 2 additional velocity (VELAD). 


RETURN 
87 CONTINUE 


Ca ea allel SECTION 7 ---- VOLUMETRIC SOURCE FOR GALA 
RETURN 
88 CONTINUE 
a ta SECTION 8 --- CONVECTION FLUXES 
RETURN 
89 CONTINUE 
C I SECTION 9 --- DIFFUSION COEFFICIENTS 
RETURN 
810 CONTINUE 
C Ww mw me mew en nen == SECTION 10 --- CONVECTION NEIGHBOURS 
RETURN 
811 CONTINUE 
C eee = SECTION 11 --- DIFFUSION NEIGHBOURS 
RETURN 
812 CONTINUE 
C Wo mmm n nnn nnn mewn nn ne SECTION 12 --- LINEARISED SOURCES 
RETURN 
S135 CONTINUE 
C Wo ww mw nnn ewe ween ne Seo LON 1S] ---=CORREG PPON GOLEPIGIUENTs 
RETURN 
814 CONTINUE 
C Wa am ewww ween nenn- SECTION 14 --- USER'S SOLVER 
RETURN 
es) CONTINUE 
C ee = ee SEC MIO Sloe --- OC ANGE SOLUTION 
RETURN 
C * bea ay aad SAUD ae senenees in Spon oe 
Cwsisvevestves Cah Aetsiots Chita latela irik tat ie ts emilee iy ioc ter isclivinciae ero a tracise ees 
C 
C--- GROUP 9. Properties of the medium (or media) 
& 
C The sections in this group are arranged sequentially in their 
C order of calling from EARTH. Thus, as can be seen from below, 
C the temperature sections (10 and 11) precede the density 
C sections (1 and 3); so, density formulae can refer to 
C temperature stores already set. 
9 ae es en 92 93, fee ie oe rae sae eas 900, ,901, 902, ee ioe 
Cvrdevers Jevedesevedze 
900 CONTINUE 
le SoG. LON = == = = - -- 
C Bera LE GRND = === —===— phase-1l temperature Index AUX(TEMP1) 
RETURN 
901 CONTINUE 
C e-em ewe ~~ Ose) ee a — — — 
C Someile 2 bE GRND<<-+=<«'=> phase-2 temperature Index AUX(TEMP2) 
RETURN 
ee CONTINUE 
BT SiC DON 2 ee a = = SS 
C ie Peer on 0 Sa e phase-1 length scale Index AUX(LEN1) 
RETURN 
o> CONTINUE 
GC FW acne eweeereeecen--- SS CPCI I) I i ie 
C oe EL2. LE. GRND--------- phase-2 length scale Index AUX(LEN2) 
RETURN 
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91 CONTINUE 


CK ween nee nee ene eee--- SECTION 1 concer ccc cr ence neennenn--- 
c For RHO1. LE.GRND--- density for phase 1 Index AUX(DEN1). 
965 CONTINUE 
RETURN 
92 CONTINUE 
CF were nee erence ecene-- SECTION 2 ceccceseccccen cence rre----- 
C For DRH1DP. LE. GRND--- D(LN(DEN))/DP for phase 1 (D1DP). 
RETURN 
93 CONPINGE 
CF ewe een cence en------ SECTION 3 cece crrccscecsessnccccsasas 
C For RHO2. LE.GRND--- density for phase 2 Index AUX(DEN2). 
RETURN 
94 CONTINUE 
c Wo ere www wee ene neeenn- SECTION 4 ence cn nnn cc cere cece eecen= 
C For DRH2DP. LE.GRND--- DC LNCDEN))/DP for phase 2 (D2DP). 
RETURN 
95 CONTINGE 
OR a SECTION 5 cecccesren cen neeceee<------ 


C For ENUT. LE.GRND--- reference turbulent kinematic viscosity. 
CALL FNS5CAUX(VIST) , AUX(CTEMP1) ,4. 886,5.432E-2,-2. 426E=-5,7. 931E=9% 
l-1. 104E-12) 
CALL FN25( Aus VIST), 1. Uma o») 
CALL FN27( AUXCVIST) , AUX( DEN1) ) 
CALL FN25( AUX(CVIST) , 20000. 0) 


RETURN 
Ze CONTINUE 
Oh Sao So Same Soe El OO Ne 
C Poe ENUL. LE. GRND--- reference laminar kinematic viscosity. 
RETURN 
97 CONTINUE 
C FF mee emwwnenwecceecae SECTION 7 S28 encccrn cc cc ces csseseecsce 
C For PRNDTL( ).LE.GRND--- laminar PRANDII nos. , ot ditiusiv ies 
RETURN 
98 CONTINUE 
C een een een eee ececece SECTION 8 cece ercrrecccccccccccesce= 
C For PHINT( ).LE.GRND--- interface value of first phase(FII1). 
RETURN 
99 CONTINUE 
Ce SECTION 9 cree ere nrc ene nn meen ren HH-=- 
G For PHINT( ). LE.GRND--- interface value of second phase(FII2) 


RETURN 
Cove ders eve dete sede ce os oe ve cove Tete de B59 FS Fe Fe Te TEE eT TS Fe To ETE TE TS EE TE BE TEE TE TET ETE ETE TE ETE TEE TSE TOTS TELE TE TESTE TE 
C 
C--- GROUP 10. Inter-phase-transfer processes and properties 
C 
10 GO TO @lGietlg2 10a lew eras 
101 CONTINUE 


Co a SECTION ne srt ioe 
C For CFIPS. LE.GRND--- inter=phase friction coeff SAUX  INTERE). 
RETURN 
102 CONTINUE 
COS. Sie a el mia SECTION 2: S39 ip eee ini eae 
C For CMDOT. EQ. GRND- inter-phase mass transfer Index AUX( INTMDT) 
RETURN 


104 


103 CONTINUE 


0 er ee mmm een eo -------- SECTION 3 crrcrerc creer ete err rence ee- 

C PoreelNi(  ). EQ. GRND=-- phasel-to-interfadee transfer 

c coet ficients -eOl) 
RETURN 

104 CONTINUE 

C HW mem eww cwn sen nesecce SECTION 4 eeeccrcern nnn ccc enc n cc cceen- 

C For CINI( ).EQ.GRND--- phase2-to-interface transfer 

C eoet fictents, (Col) 
RETURN 

Sos) aria tr oiroke Sri ice bir bie ark tee Se ae eh eae ee areal a keh bi eee eka ck Dr ae ak arbi ae sear bak kaka 


C 
Eee- GROUP 11. Initialization of variable or porosity fields 
C 
11 CONTINUE 
RETURN 
Cretedeteve seve sede tose se te Ve TEPC T ER TE TE TE TE TE TE CTE TEETER TENE TE TE TENE FETE ETE TE SE TELE TENE TER EE TEE TERETE TENE TEE TEC 
C 
C--- GROUP 12. Convection and diffusion adjustments 
C 
12 CONTINUE 
RETURN 


(epee 6. § 6 6 CaCI PoC PURPA Teak Joel Teeal Sect TO aC Herma goat ae one te Jo ceel Jena ral eal Wel pyre eT et SU ead ee one Joes epee Yoana, Pe Pee teatetis octal ts alto estciects atin ta alan latent ce twa” 
Cwseuveves LL Le ee Lr) Sele ark eke) rl ek rd de Le ed ek ek ee ek ee rd ie ed ee ee) Se ek ed ee | 


% 
sip kar hors irk irk rs ers er) ) Ord ek) er | err ees 


C--- GROUP 13. Boundary conditions and special sources 
C 
13 CONTINUE 
eS Omwemers® ,131,1325133,1345135, 136,137 , 135,139 ,1310, 
immcminwee de? 1919 13 16.1315 dio 37 9318 1919 ,1320,1321),1sC 
130 CONTINUE 


Coen n en nce nee ceee SECTION 1 -eerccccccee- coefficient = GRND 
131 CONTINUE 
Corer ren ene meee eee SECTION 2 erccccer----- coefficient = GRND1 
RETURN 
ime CONTINUE 
C-------- eee eee SECON ee = coefficient = GRND2 
RETURN 
133 CONTINUE 
C------------------- SECTION 4 --crrerccceee coefficient = GRND3 
Rel URN 
134 CONTINUE 
Cone --------- eee -- Sed l¢Na. > ---.--------— coefficient = GRND4 
RETURN 
135 CONTINUE 
ee = = == SU Oil ONO syeaie ea a =e oe coefficient = GRND5 
RETURN 
136 CONTINUE 
Oat SEURION: 7 ----<--------- coefficient = GRND6 
RETURN 
137 CONTINUE 
C------------------- SECTION 8 -rrrrererrnc- coefficient = GRND/ 
RETURN 
138 CONTINUE 
Coe --- oe --------- SECTION 9 e------------ coefficient = GRND8& 
RETURN 
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139 CONTINUE 


EE SECTION 10 ------------- coefficient = GRND9I 
RETURN 
1310 CONTINUE 
Oe SECTION 11 ------------- coefficient = GRND10 
RETURN 
1311 CONTINUE 
ee eae SECTION 12 -------------------- value = GRND 
1351 460N TENCE 
RETURN 
1312 CONTINUE 
ee ee ee ee = SECTION 13 -------------------- value = GRND1 
RETURN 
1313 CONTINGE 
Cos -ceeeeee eee ae SECTION 14 -------------------- value = GRND2 
RETURN 
1314 CONTINUE 
en ee SECTION 15 -------------------- value = GRND3 
RETURN 
1315 CONTINUE 
Os < na Becta eels o SECTION 16 -------------------- value = GRND4 
RETURN 
1316 CONTINUE 
C+ - << 6822225 eee eee SECTION 17 -------------------- value = GRND5 
RETURN 
1317 7COh TINUE 
ee eee ee SECTION 18 -------------------- value = GRND6 
RETURN 
1318 CONTINUE 
ee eee ae SECTION 19 -------------------- value = GRND?7 
RETURN 
1319 CONTINUE 
On=- (22ers SECTION 20 -------------------- value = GRND8 
RETURN 
1320 CONTINUE 
eee ne aa SECTION 21 -------------------- value = GRND9 
RETURN 
1321 CONTINCE 
ee Se eae SECTION 22 -------------------- value = GRND1O 
RETURN 
Credevedeseseds reve leds e Te sete Ne se ted Ve Fe eee TONE TE ETE TOT Te ENE TENE ENE TE TED E TE TENE TONE TE VEN TE TE TENS TET TENE TEN TC TE TCG 
C 
C--- GROUP 14. Downstream pressure for PARAB=. TRUE. 
C 
14 CONTINUE 
RETURN 
Cvevevevedevedese cede de Bevel Tee Tete Te Tee TS TS Te ETE TE TOTS TSS ISTE TS ESTE TS TS TS TSEC TS TTC TEE TE TTR EERIE RACER 
C 
C--- GROUP 15. Termination of sweeps 
C 


15 CONTINUE 
C * Make changes for this group only in group 19. 
RETURN 


mlaale ale atealan!ecianta wlentaalon’sntsnlon'ontoa'salen'saleetaatsatantaa'saiantnnlantaalasntan'sn sntsaisa'saisuisatea'snioa's ale alentsntan'entanteulentan! sleslesion’ elaalsalsatou'an'te s'e 
EN ER ER ER ED ED EN ES EL ER ER ER VEN FR ER EN FL ER ED ED FL ECR FR FRED GR PL ER FL GR GL ER FR ER FB EL GR ER EL CHR EL ER ER ER EN CL FT Hd FR FD 7 #xerv en riverty eu ed rierts ev¢ 


C 
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w--- GROUP 16. Termination of iterations 
C 
16 CONTINUE 
C * Make changes for this group only in group 19. 
RETURN 
Cr vedessseds sede seve sete de ede ve se te ede Bese Neve ve Seve Te Fete Fe Fe Tete Te Ve Te Ve TET Tete Ne Tes Tee Te Tes ede secede setetesede sete 
C 
C--- GROUP 17. Under-relaxation devices 
C 
17 CONTINUE 
C * Make changes for this group only in group 19. 
RETURN 
Credededsdedevs ve dese dete te seds tote de Teas Fete Tete ve Fe Fe Pete TeV Te Te TE ETE V ETE TE TENE VEE TE TEV VE TEV VED TETE NTE TeV TTT TE TeV 
C 
C--- GROUP 18. Limits on variables or increments to them 
C 
18 CONTINUE 
Gee Make changes for this group only in group 19. 
Re TURN 
Creve dedede te tere de Te ve Ne te Fe Fe e8 Vere Te Fe TENE TE TET TENE TENE TE TE Te TE TE TE TEE TS ETE FE TENE TE TENE TEE TE TE FETE TE GEE PETE TE TE TENE TE 
C 
G=== GROUP 19. Special calls to GROUND from EARTH 
C 
Meee TO (191,192,193,194,195, 196,197,198), 1SC 
m1 CONTINUE 


2 i ee SECTION 1 ---- START OF TIME STEP. 
RETURN 
192 CONTINUE 
C  * --+---------------- OECD ION ere oo ahd Or oWEED: 
RETURN 
mS CONTINUE 
Cl SECTION 3 ---- START OF IZ SLAB. 
RETURN 
194 CONTINUE 
 » 4 ae SECTION 4 ---- START OF ITERATION. 
RETURN 
is CONTINUE 
ke SECTION 5 ---- FINISH OF ITERATION. 
RETURN 
Ho) CONTINUE 
Wo weer ween eee eee ee ee ee SECTION 6 cS FINISH OF IZ SLAB. 


OPEN (UNIT = 10, FILE = 'DATA’) 

Se eGETYx( ¥1.GV1 ,MY , MX) 

CALL GETYX(W1,GW1,MY,MX) 

Ait GEryX( AUX(TEMP1) ,GIMP1,MY,MX) 
CALL GETYX(AUX(DEN1) ,GRHO1,MY,MX) 

CALL GETYX(H1,GH1,MY,MX) 

CALL GETYX( AUX(VISL) ,GENUL,MY ,MX) 
Ho1000 1=1.MY 

GSONIC = (GAMA*RAIR**GTIMP1(1,1))**. 5 
GVEL = (GV1(1,1)**2 + GW1(1,1)**2)**. 5 
iWaeci see EG. 10) . AND. (1Y -EQ: 1) / AND: 
1(ISWEEP .EQ. LSWEEP))THEN 

WRITE(10,*) 'GV1 =', GV1(1,1) 
Weigeaiee)\' GIMP] = , GIMP1(1;1) 
WRITE(10,**) 'GW1 =', GW1(1,1) 


Cy Cara ae a) a 2G 6) GG 
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END IF 

GMACH(I,1) = GVEL/GSONIC 

CALL SETYX(C1,GMACH ,MY,MX) 

CALL GETONE( AUX(TEMP1),TINF ,MY, 1) 

CALL GETONE(C1,MINF ,MY, 1) 

CALL GETONE(W1,WINF ,MY, 1) 

CALL GETONE( AUX(DEN1) ,RHOI ,MY, 1) 

GTR = TINF*(1 + (RFAC*(GAMA - 1. ))**(MINF?*2)/2) 
DO 1010 I=1,MY 


Cc 
a 
OQ 
S 
© 


GHE( 1.1) = GO( 11) (Gira Wal 
GST(I,1) = GHC(I,1)/(RHOI*WINF*CP) 
GTAU(I,1) = GRHO1(1,1)*GENUL(I,1)*GW1(1,1)/2. 5E-4 
1010 GCF(I,1) = 2*GTAU(I,1)/CRHOI*WINF**2) 
CALL SETYX(C2,GHC,MY,MX) 
CALL. SETYR( G3, GS ayo 
CALL SETYX(C4,GCF,MY,MX) 
CLOSE (UNIT = 10) 
RETURN 
197 CONTINUE 
Co ----------- +--+ eee SECTION 7 ---- FINISH OF SWEEP. 
RETURN 
198 CONTINUE 
Cw -ee---------------- SECTION 8 ---- FINISH OF TIME STEP. 
RETURN 


a 2e* 26 26 Sous a's Ae alfa ale as ao ate ae ats wants afsals ats aloats ale ale o's ala ale als ae ale ate els ale alo ats ats a's ale we ale ate ate al onts ale ee ale ele awe ale Se ale ele ate wle 
ard ark ard er ard ard ar ed er rd ed ek rk ed ed ee de a ek ek rd ek ek er ek es ek kk ee ed ed a ed ed ed a er re ek a er ek er 


Ca Gee Cao CGa-C C2 CoC) CCC ca 


C--= GROUP 20. Prelimanary epmiae one 
G 
20 CONTINUE 
RETURN 

Crevevevedsvedeceds doves 
GC 

C==- GROUP 21. Printeout of Wwarwables 
C 


we alanla alone al salontsnlonts alon'salon'cnlontontwalonlontononion!oalon' suleoniontiontantontoniontontontoatontonlon'salsaton calentontontss'se'oalsatonte 
CN ON ON ARAN GN GLEN GL CC GL BRL ERED EN SL GC OR GO GLEN BBN EL EVOL ERED EV GRE TLON FC GV GL ET IN GEL GLEN EVEL SGN ED ON GN GS 


21 CONTINGE 
C * Make changes for this group only in group 19. 
RETURN 
Okikieiete thik tk reat ete 
C 
G-=- GROUP®22.) Spoue-val WweNprime come 
22 CONTINUE 
C * Make changes for this group only in group 19. 
RETURN 
Crsvedededsdodeveds seve 
C 
C--- GROUP 23: Field printout andgpieemcenc.c! 
23 CONMINGE 
RETURN 
Orikteeeieibieisi¢iete a teihitikibihit tite titieieie site ee eet hie ee te ee ek eee ee ee ee ee ee ees 
C 
C--- GROUP 24. Dumps for restarts 
C 


meee teste tetesesetededededesesedededededesedededvededsededcdedesetevetededededededesevesideve 


Soulontectontioctonionteclonlont ents alentsnlonts clsniscten'on'selentonlontoniacicctonts ale wtsnlsatoals lentes alsctsclentontoniostoctssloulontonlaataalen'sate 
GR ARAL ON AN GR EL GLEN GL ER ALOR AD FL AL GD GL AL IL GL AL OL OL EN GN GL IN ARGS ERAN ENGR ENGLER AD ON GN ARAL IL ER AS ON GR ER OL ON ON ON 


24 CONTINUE 
RETURN 
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GQ(1I,1)=(GRHO1(1,1)*GENULCI,1)/. 7)*((GH1( 1,1) -CP=TWAn)y 23 poem 


END 
i v 
// 


B. COMPGRD 

//COMPGRD JOB (2529,9999), 'NIEBURG' ,CLASS=B ,REGION=2000K 
//*MAIN LINES=(99) 

fei EXEC FORTVCL, 

yy PARM. LKED=( ' LET ,NCAL,SIZE=(5000K, 64K) ,NOXREF' ) 
//FORT.SYSIN DD DSN=MSS. 81541. PHOENICS. SRC(GROUND) , DISP=SHR 
//FORT. SYSLIB DD DSN=MSS. PHOENICS. COMMONS. FORT, DISP=SHR 
//LKED.SYSLIB DD 

ii DD 

Wi DD 

//LKED. SYSLMOD DD DSN=MSS. $1541. PHOENICS. LOAD( GROUND) , DISP=SHR 


C. LINKEAR 

//LINKEAR JOB (2529,9999),'*CHAM' ,CLASS=B ,REGION=7000K 

//*MAIN LINES=(99) 

//LINK EXEC PGM=IEWL, PARM=' XREF , LET, LIST,SIZE=( 2000K, 256K) , AMODE=24 ' 

V/eYoPRINT DD SYSOUT=* 

P7eyoUTl DD UNIT=SYSDA, SPACE=(CYL,(1,1)) 

Py =siSoLIB DD DSN=MSS. $1541. PHOENICS. LOAD, DISP=SHR ,MSVGP=PHOENICS 
DD DSN=MSS. PHOENICS. EARTHL. LIB, DISP=SHR,MSVGP=PHOENICS 
DD DSN=MSS. PHOENICS. RESLIB. LIB, DISP=SHR,MSVGP=PHOENICS 
DD DSN=MSS. PHOENICS. SATLIB. LIB, DISP=SHR,MSVGP=PHOENICS 

DSN=SYS1. VLNKMLIB, DISP=SHR, UNI T=3350 

DD DSN=SYS1. VFORTLIB , DISP=SHR, UNIT=3350 

SYSLMOD DD DSN=MSS. 581541. PHOENICS. LOAD, DISP=SHR,MSVGP=PHOENICS 

SYSLIN DD * 

INCLUDE SYSLIB( GROUND) 

INCLUDE SYSLIB(GREA1) 

INCLUDE SYSLIBCMACDEP84 ) 

ENTRY MAIN 

NAME EARTHP(R) 


// 


SS 
ee oe 
i) 

SS 
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APPENDIX E. PROGRANIS REQUIRED FOR COMPACTING FILES 


A. PHOBNX0!1 LOADLIB 
//RATNER JOB (2529,9999),' PHOENIX LOAD LIB’ ,CLASS=A 
//* STEP 1 RENAME CURRENT PHOENIX LOAD LIBRARY TO TEMPORARY LIBRARY 
//STEP1 EXEC PGM=IDCAMS 
//SYSPRINT DD SYSOUT=A 
//SYSIN@DD * 
ALTER “MSS.$1541.PHOENICS. LOAD - 
NEWNAME(MSS. $1541. PHOENICS. LOADTEMP) 
//* STEP 2 AULOCATE SPACE FOR NEW PHOENIX LOAD LIBRARY 
//STEP2 EXEC PGM=IEFBR14 
//DD1 DD UNIT=SYSDA,DISP=(NEW, CATLG, DELETE), 


wy) DCE=(RECEM=U ,PLKSmzE—13050)), 
1) SPACE=(13030,(400,15,7)), 

iy DSN=MSS. $1541. PHOENICS. LOAD 
/ os 


B. PHOEN X02 LOADLEE 
a 
//RATNER JOB (2529,9999),'PHOENIX LOAD LIB' ,CLASS=A 
//* STEP 3 COPY RENAMED OLD LOAD LIBRARY TO NEW LOAD LIBRARY 
//STEP3 EXEC PGM=IEBCOPY 
//SYSPRINT DD SYSOUT=A 
//IN DD DSN=MSS. $1541. PHOENIGS. LOADTEMP ,DISP=(OLD, KEEe) 
//OUT DD DSN=MSS.S1541. PHOENICS. LOAD, DISP=(OLD, KEEP) 
//SYSIN DD * 
COPY EB DD=IN , OUTDD=OUT 
/* 
// 


C. PHOEMN0S LOADLIB 


//RATNER JOB (2529,9999),' PHOENIX LOAD LIB' ,CLASS=A 
//* STEP 4 DELETE REWAMED OLD PHOEBE SI Ss EOay wiih Ran 
//STEP4 EXEC PGM=IDCAMS 
// Sisee ea Dies, sOUT—A 
/ {Sie 
DELETE 2155. 5154). PHORNIGse, LOsDoEE NONVSAM SCRATCH 
7 ry 
// 
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D. PHOENNO! SOURCELB 


Pe eR JOB aoe 9999) peda aca ees LIB | eee i 


ale s'es‘e sles ‘os! tewe eset esis nfo sectestes" evesies'es' es'es's ate a's ae ss = eve wee! efenes 
ay ray Boos ev even exv eves 


= 


ee STEP 1 RENAME, CURRENT PHOENIX ‘SOURCE LIBRARY TO TEMPORARY LIBRARY 


UNS ae aN aN ON EN ES a inl Wane Siri =. as ane evesss Ay Phy Pats 


— 


VI STEP] EXEC PGN=IDCAMS: 

//SYSPRINT DD SYSOUT=A 

//SYSIN DD * 
ALTER MSS. $1541. PHOENICS. SRC - 
_NEWNANE( NSS. $1541. PHOENICS. SRCTEMP) 


I s'es'e ses! s'es's eyes le es oe! sles =" ewe eslesies's eves's Je s'e tevle s'e Je Je s"¢ estesle s'e se s'e s'< eve s'e Je slevJe lesies'e ate s'e ou! ss este wJe ate Je s' ese Je Jee ~ 4 ewe Je wledantsctswasntawaselsulsn'sese'aa's 
es es ky ER ERM ETAL IR ER ERM GL Gd IV Ev ErN ED 


//* STEP 2 “ALLOCATE SPACE FOR NEW PHOENIX SOURCE LIBRARY 

//* * wevetededs evs See teledede eset Cece te Cece Te Te ee Te TEP ONG SOLOS ONE TE TET VOTE CE TET TEDODE TEE TESTE TOPE GETE 
/ /STEP2 EXEC PGNSIEFBR14 

//DD1 DD UNIT=SYSDA,DISP=(NEW,CATLG, DELETE) , 

fo DCB=( RECFN=FB , LRECL=80 , BLKSIZE=800), 

SPACE=( 23440 ,(200,30,7)), 

DSN=MSS. $1541. PHOENICS. SRC 


SOS SS 
— oo 


FE. PHOENX02 SOURCELB 


a. Pie Rees 9999) | Seas seenciraeas LIB ees is 

fe 2S leila iv ig Gira k cs orig ks Ties riers ntsc oer sila arnchars 
/* 
es 


Yaw en's non’, a atae an’ wee! anton! ae’, tam tevlesesles als nlenlaatan's 
Shiner crn viele i wriwiy v% eveve 


| oe EXEC PGM=IEBCOPY | 
V7 o> 1oPRINT DD SYSOUT=A 
oie DD DSN=NSS.$1541. PHOENICS. SRCTEMP, DISP=( OLD , KEEP) 
[7 oP DD DSN=NSS. $1541. PHOENICS. SRC ,DISP=(OLD, KEEP) 
PyersiN DD * 
COPY INDD=IN,OUTDD=OUT 
/* 
// 


F. PHOENNO3 SOURCELB 


//RATNER JOB (2529,9999),' PHOENIX SOURCE LIB' ,CLASS=A 
V7 STEP 4 DELETE RENAMED OLD PHOENIX SOURCE LIBRARY 
//STEP& EXEC PGM=IDCAMS 
//SYSPRINT DD SYSOUT=A 
//SYSIN DD * 
Merete MSS.81541. PHOENICS. SRCTEMP NONVSAM SCRATCH 
/* 
// 
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APPENDIX F. PROGRAMI “FORNIO9” 


"FORMO9" is required to format the PHOENICS output file DFO9 to @ files 


can be read by the graphics program “PITOTON’. 
//FORMO9 JOB (2529,9999), ‘CHAM’ ,CLASS=B 
//*MAIN LINES=(20) 

// EXEC FORTVCLG 

// FORTS Youn. Dir: 


C PROGRAM FORMO9 
C 
C Increase PARAMETER settings if more variables or a larger grid 
C is required. 
C 
PARAMETER( NPHI P=80 ) 
C NB NXYMAK must be >= NZ and NX*NY 


PARAMETER( NX YMAK=20000) 

CHARACTER*4 NAMECNPHIP) ,MESS( 10) 
CHARACTER=15 FILIN, FILOUT 
CHARACTER*1 ANS 

CHARACTER*10 ACCESS 

LOGICAL LBUF(10),STORECNPHIP) , DEBUG 
DIMENSION RBUFCNXYMAX) , IBUF( 10) 


100 FORMAT(1X,71I10) 
101 FORMAT(6(1PE13. 6)) 
102 FORMAT(1X,79L1) 
103 FORMAT(1X,19A4) 


on®) 


Reset DEBUG=. TRUE. to get error messages.. 
DEBUG=. TRUE. 


I SB 


set. Logical sUnibs. 2. 
LUIN=9 
LUOUT=10 
LUBUG=6 


Set File Names... 

Name of unformatted restart file 
WRITE(6,*) ' RESTART FILENAME? ' 
READ(5,'(A)') FILIN 

Name of formatted restart file 
WRITE(6,*) ' FORMATTED FILENAME?’ 
READS. CA) aeour 


Direct access or sequential? 
WRITE(6,*) ' DIRECT ACCESS (PHIDA) FILE? Y/N’ 
READ(5,'(A1)') ANS 
IF(ANS. EQ. 'Y'. OR. ANS. EQ. 'Y') THEN 
Version 24 cera t 
ACCESS=' DIRECT’ 


ey ec een) Court Cea GC) Cac 3 GC) 


C 
C 


LENREC=256 
ELSE 


Ssversion 1.3 


CF CG Garey Cig O) 


>) 


C 


ACCESS=' SEQUENTIAL ' 

LENREC=100 

ENDIF 

WRITE(6,*) ‘ DEFAULT LENREC? (100 V1.3, 256 V1.4) - Y/N' 
READ(5,'(A)') ANS 

IF(ANS. NE. 'Y'. AND. ANS. NE. 'Y') THEN 

WRITE(6,*) ' LENREC? (WORDS NOT BYTES)' 

READ(5,'(1)') LENREC 

ENDIF 


NBYTES is the number of bytes per single precision word. 


NBYTES=4 


Record length in bytes. 


LREC=NBYTES*LENREC 


SeeOrEN input file: 


C 


OPEN( UNIT=9 , FORM=' UNFORMATTED' ,STATUS='OLD' , 


1 ACCESS=ACCESS , IOSTAT=IERR , ERR=90 ) 


IF(IERR. NE.0) GO TO 90 
IF( ACCESS. EQ. 'DIRECT') THEN 
READ( LUIN,REC=1) (MESS(I),I=1,10) 
ELSE 
READCLUINe( MESS(1),1=1,10) 
ENDIF 
IF(DEBUG) WRITE(LUBUG,*) MESS 
MCACEESS 60, DIRECT’) THEN 
READ( LUIN,REC=2) (LBUF(1),I=1,4) 
ELSE 
READ(LUIN) (LBUF(1),I=1,4) 
ENDIF 
IF(DEBUG) WRITE(LUBUG,*) (LBUF(I),1I=1,4) 
IF( ACCESS. EQ. ‘DIRECT’ ) THEN 
READ  LUIN. REG=3) NX,NY,NZ,NPHI,(IBUF(1),1=1,6),LENREC 
ELSE 
READ(LUIN) NX,NY,NZ,NPHI,( IBUF(1I) ,I=1,6),LENREC 
ENDIF 
NXNY=NX*NY 
IF(NXNY. GT. NXYMAX) GOTO 92 
IF(NPHI. GT. NPHIP) GO TO 93 
IF(DEBUG) WRITE(LUBUG,**) NX,NY,NZ 


Peper BNe output file: FORMATTED, SEQUENTIAL and of FIXED RECORD LENGTH. 


OPEN( UNIT=10, FORM=' FORMATTED’ , STATUS='OLD' , 


1 ACCESS='SEQUENTIAL' , IOSTAT=IERR,ERR=91) 


POOIe kh Rane) GO) TO 91 
Wiis CiUOUT, 103) (NMESS(1),1=1,10) 
Mette CLUOUT 102) CLBUF(I),I=1,4) 
WRITECLUOUT,100) NX,NY,NZ,NPHI,(CIBUFCI),1I=1,6) , LENREC 
IF( ACCESS. EQ. ‘DIRECT’ ) THEN 
READ( LUIN,REC=4) RINNER 
ELSE 
READCLUIN) RINNER 
ENDIF 
WRITE( LUOUT,101) RINNER 


LIS 


IF( ACCESS. EQ. ‘DIRECT’ ) THEN 
READ( LUIN, REC=5) (NAMEC I) ,I=1,NPHI) 


Eeob 
READC(CLUIN) (NAMECI) ,1=1,NPHI) 
ENDIF 
C Set C36... names for variables 51 to NPHI for version 1.4, touarema 


C control characters etc| (Unless already stored) 
IF( ACCESS. EQ. 'DIRECT'. AND. NPHI. GT. 50) THEN 
READ( LUIN, REC=10) (STORE(1),I=1,NPHI) 
DO 6 I=51,NPHI 
J=1-15 

IF(. NOT. STORE(1)) THEN 

IF(I. LT. 100) THEN 

WRITECNAMEGI), © © a e2qen se” @) ones) 

ELSE 

WRIDECNAME@I =. @ (Gl oeite aes] 

ENDIF 

ENDIF 

6 CONTINUE 
ENDIF 
WRITE( LUOUT,103) (NAME(I),I=1,NPHI) 
IF(DEBUG) WRITE(LUBUG,*) (NAME(I),I=1,NPHI) 
IF( ACCESS. EQ. 'DIRECT') THEN 
READ( LUIN, REC=6) (RBUF(I),1I=1,NX) 
ELSE 
READ(LUIN) (RBUF(I),I=1,NX) 
ENDIF 
WRITE( LUOUT, 101) (RBUF(I),I=1,NX) 
IF( ACCESS. EQ. 'DIRECT') THEN 
READ( LUIN, REC=7) (RBUF(I1),I=1,NY) 
ELSE 
READ(LUIN) (RBUF(I),I=1,NY) 
ENDIF 
WRITE( LUOUT,101) (RBUF(1I),I=1,NY) 
IF( ACCESS. EQ. 'DIRECT' ) THEN 
READ( LUIN,REC=8) (RBUF(1I),I=1,NZ) 
ELSE 
READ(LUIN) (RBUF(I),I=1,NZ) 
ENDIF 
WRITE( LUOUT,101) CRBUF(1),1I=1,NZ) 
IF( ACCESS. EQ. DIRECT’) THEN 
READ( LUIN, REC=9) (RBUF(1),1I=1,NZ) 
ELSE 
READ(LUIN) (RBUF(I),I=1,NZ) 
ENDIF 
WRITE( LUOUT,101) (RBUFCI),I=1-NZ) 
IF( ACCESS. EQ. 'DIRECT') THEN 
READ( LUIN, REC=10) (STORE(I),I=1,NPHI) 
ELSE 
READ(LUIN) (STORE(1I),1I=1,NPHI) 
ENDIF 
WRITE( LUOUT,102) (STORE(1),I=1,NPHI) 
IF(DEBUG) WRITE( LUBUG,102) (STORE(1I),1I=1,NPHI) 

NVAR=0 
DO 5 =i Neen 
IF(STORE(I)) NVAR=NVAR+1 


li4 


5 CONTINUE 

TREC=10 

DO 10 K=1,NZ 

DO 20 I=1,NVAR 

IF( ACCESS. EQ. DIRECT’) THEN 

IREC=IREC+1 

NVRECHNXNY / LENREC 

IF(MOD(NXNY , LENREC). NE. 0) NVREC=NVREC+1 

DO 30 IVREC=1,NVREC 

or=( 1VREG-1)*LENREC + 1 

JL=IVREC*LENREC 

JL=MINOCJL,NANY) 

READ( LUIN,REC=IREC) (RBUF(J),J=JF,JL) 
30 CONTINUE 

GSE 

READ(CLUIN) (RBUF( J) ,J=1,NXNY) 

ENDIF 

WRITE( LUOUT ,101) (RBUF(J) , J=1,NXNY) 
20 CONTINUE 
10 CONTINUE 

Pe DEBUG) Wi CLUBUG 1900 ) 
1900 FORMAT(1X,' NORMAL STOP IN PROGRAM’ ) 


GO TO 999 
Seercor Crapping... 
90 IF(DEBUG) WRITE( LUBUG, 1901) 
1901 FORMAT(1X,' ERROR IN OPENING INPUT FILE') 
GO TO 999 
91 IF(DEBUG) WRITE( LUBUG, 1911) 
1911 FORMAT(1X,'ERROR IN OPENING OUTPUT FILE' ) 
GO TO 999 
92 IF(CNXNY. GT. NXYMAX) WRITE( LUBUG, 1921) 


DC NZ5GT.NXYMAX): WRITE( LUBUG. 1923) 
1921 FORMAT(1X,'SETTING OF NX*NY TOO LARGE; RESET IN FORMO9 AND ' 
1,/,1X, 'RE-CREATE THE LOAD MODULE’ ) 
1923 FORMAT(1X,'SETTING OF NZ TOO LARGE; RESET IN FORMO9 AND ' 
1,/,1X,'RE-CREATE THE LOAD MODULE’) 
GO TO 999 
93 WRITE( LUBUG, 193) 
193 FORMAT(1X,'SETTING OF NPHI TOO LARGE; RESET IN FORMO9 AND ' 
1,/,1X,'RE-CREATE THE LOAD MODULE’ ) 
999 CONTINUE 
MmEGhOSE files ... 
SHOSEC LUIN ,STATUS="KEEP' ) 
CLOSE( LUOUT, STATUS=' KEEP' ) 
STOP 
END 
/GO.FTO9F001 DD DSN=MSS. $1541. PHOENICS. AUG90, DISP=( OLD) 
/GO. FT10F001 DD DSN=MSS. $1541. PHOENICS. DFO09C, DISP=( OLD) 


”~ 
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